Cold Dark Clouds 



1 



o 



> 
in 



X 



COLD DARK CLOUDS: The Initial 
Conditions for Star Formation 



Edwin A. Bergin 



, Department of Astronomy, University of Michigan, 500 Church St. Ann Arbor, 



MI, 48109, USA; email ebergin@umich.edu 



i3 : Mario Tafalla 

Observatorio Astronomico Nacional, Alfonso XII 3, E-28014 Madrid, Spain; 
1/-^ . email: m.tafalla@oan.es 

(N 



Key Words interstellar medium, interstellar molecules, molecular clouds, pre- 



O ' stellar cores, star formation 



Abstract Cold dark clouds are nearby members of the densest and coldest phase in the 
galactic interstellar medium, and represent the most accessible sites where stars like our Sun are 
currently being born. In this review we discuss recent progress in their study, including the newly 
discovered infrared dark clouds that are likely precursors to stellar clusters. At large scales, dark 
' clouds present filamentary mass distributions with motions dominated by supersonic turbulence. 

, At small, sub-parsec scales, a population of subsonic starless cores provides a unique glimpse 

fT^ ' of the conditions prior to stellar birth. Recent studies of starless cores reveal a combination 

I of simple physical properties together with a complex chemical structure dominated by the 

^— «j . freeze-out of molecules onto cold dust grains. Elucidating this combined structure is both an 

' observational and theoretical challenge whose solution will bring us closer to understanding how 

, molecular gas condenses to form stars. 
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1 INTRODUCTION: HOLES IN THE HEAVENS 

"Hier ist wahrhaftig ein Loch im Himmel!" or "Here is truly a hole in the heav- 
ens!", William Herschel was heard by his sister Caroline to exclaim in 1784 when 
he pointed his telescope towards the con stellation Scorpius and viewed the dark 
nebulae in Qp hiuchus (iHoughtonI 19421 ). Herschel reported his discovery the 



following year fascy [izd ) , an ^lTr the next century astronomers debated 



whether these dark objects where truly voids, left by the stars as theorized by 
Herschel, or were perhaps nebulous dark regions observed projected upon a dense 
and bright stellar background. Our views began to change when Edward Emerson 
Barnard, who was motivated by Herschel's discovery, published the firs t modern 



and s ystematic photographic survey of the "Dark Markings of the Sky" (jBarnard 



mi). Barnard argued that his deep photographs provided increasing evidence 



that many of these dark areas were "obscuring bodies nearer to us than the 
distant stars." 

Such regions might remain just obscured curiosities were it not for the clear 
association with star formation which began to be recognized in the mid-poi nt of 
the t wentieth century. The first person to do so was Bart J. Bok in 1946 



19481 ) when he claimed that the dark nebula, in particular compact and nearly 



round ones that now bear the name "Bok Globules" , are the sites of stellar birth 
It was the advent of modern infrared and millimeter-wave technology in the latter 
half of the twentieth century that cemented the relation between Dark Clouds 
and the formation of stars and planetary systems. 

The discovery of molecules in space further revealed that D a.rk Clouds are 



made of molecular mate rial with H2 as the dominant constituent (IWeinreb et al 



196^ : IWilson et al]ll97nl ^ . Since this finding, the terms Molecular Cloud and Dark 



Cloud are often used interchangeably, as they refer to the two main characteristics 
of the clouds: their molecular composition and their opaque optical appearance. 
Molecular clouds are dark not because of their hydrogen molecules, but because 
of a population of tiny solids ("dust grains") that absorb the optical starlight 
and lead to high visual extinctions {Ay > l™)- Such a dimming of the starlight 
reduces the heating effects from external radiation and results in temperatures a 
few degrees above the 2.7 K cosmic background (T ~ 10 K). 

For the purpose of this review we will use the term Cold Dark Cloud to refer to 
those molecular clouds that are close enough (< 500 pc) to be seen in silhouette 
against the background galactic starlight. These clouds (e.g. Taurus, Perseus, 
Ophiuchus, Lupus, ...) are observed to be forming low mass stars either in iso- 
lation or in small compact groups and have gas temperatures ~ 10 — 20 K. This 
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stands in contrast to diffuse clouds in the Milky Way, which contain molecular 
material but are not forming stars and are not optically opaque with Ay ^ I*". 
At the other end of the spectrum of galactic star formation lie the more distant 
and more massive Giant Molecular Clouds that form rich stellar clusters and 
contain embedded massive stars that heat the surrounding gas to temperatures 
> 20 K. 

Over the past decade the development of sensitive continuum and heterodyne 
detectors at millimeter and sub-millimeter wavelengths has led to clear and im- 
portant advances in our knowledge of the physics and chemistry of these Cold 
Dark Clouds. Particular progress has been made in the study of regions prior 
to onset of star formation, which is the primary focus of this review. We will 
first discuss the gains in our knowledge starting from the larger tens to hundreds 
of parsec scale of molecular clouds. We will follow with a detailed discussion of 
the physics and chemistry of the much smaller, ~0.1 pc, scale of molecular cloud 
cores, where recent research has clarified a number of outstanding issues. We 
conclude with a discussion of a newly discovered class of cold {Tgas 10 — 20 K) 
dark clouds: infrared dark clouds that represent a new population of dense molec- 
ular regions and are believed to be the precursors to stellar clusters and massive 
stars. 

Space constraints prohibit this review from complete coverage of the forma- 
tion of clouds, stars, and planets. For m ore details we refer the reader to recent 
reviews of star formation (Larson 20031 ) (see also the review chapter by McKee 
and Qstriker, this volume) , the role of turbulence (IMac Low &: KlessenI l2004l : 

Elmegreen &: Scalol l2004l ) . pr obes of physical conditions (lEyaiis ligggl). gas/dust 

chernistry and star forma tion ( van Dishoeck and Blake 19981 : Ehrenfreund and Charnlev 
200C : Ivan Dishoecklliooi ). and the entire Protostars and Planets V volume (iReipurth et al 
20071 ). especially the papers by Ceccarelli et al., Di Francesco et al., and Ward- 
Thompson et al. that contain complementary reviews of this topic and include 
references that for reasons of limited space could not be mentioned here. Re- 
cent accounts of how the new observations with the Spitzer Space Telescope 
are advancing our k nowl edge of star forma tion in dark clouds can be found in 
Werner et aiT(|200()l ^ and lAllen et~^ (|2007l ). 

While recent progress clearly warrants this review, the year 2007 presents an 
auspicious occasion for an overview of dark clouds. It is 150 years since the birth 
of E. E. Barnard in Nashville in December 1857, 101 years since the birth of B. 
J. Bok in the Netherlands, and 80 years since Barnard's Photog raphic Atlas o f 
Selected Regions of the Milky Way was published posthumously f|Barnard 19271 ). 
It is also the centennial of Barnard's beautiful exposure of Taurus- Auriga, which 
is shown in Fig. 1. The hundred years that separate us from Barnard's picture 
of Taurus have brought us a wealth of knowledge on the physics and chemistry of 
these dark patches. As we will see in this review, however, dark clouds remain in 
some aspects as mysterious and fascinating as when Barnard pointed his camera 
to them. 



2 LARGE SCALE STRUCTURE OF COLD DARK CLOUDS 
2.1 Recent Studies 



A number of technical developments over the last decade have rapidly increased 
our ability to study dark clouds. Large-scale digital images in the optical, in- 
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frared, and radio regimes have become available in a routine or quasi-routine 
basis, and their analysis allows studying degree-sized portions of the sky at ar- 
cminute or better resolution. In the following sections we review some of these 
developments ordered by wavelength, and we briefly discuss the analysis tech- 
niques used to exploit this new wealth of data. Thanks to this recent progress, 
the study of cloud structure is currently enjoying a new golden age. 
2.1.1 OPTICAL AND IR DATA The Digitized Sky Su rvey (DSS) a nd the 
USNO Precision Measuring Machine project (USNO-PMM; iMonetl Il996l ) have 
provided multi-band, all-sky images at optical wavelengths with accurate cali- 
bration and registration. Using different techniques, these data can b e proce ssed 
automatically to create large-scale images of dark clouds. ICambresvl (| 19991 ) has 
used USNO-PMM data together with a star-counting method to produce ex- 
tinction (Ay) maps of up to ~ 250 square degrees for 24 dark clouds including 
Taurus, Perseus, and Ophiuchus. He has used an automatic adaptive grid with a 
fixed number of stars per cell to avoid empty pixels and to achieve a maximum 
extinction of around Ay ~ 7. This adaptive cell method, however, degrades the 
angular resolution of the images from about 1 arcminute in the outer parts of the 
clouds to about 10 arcminutes in the most opaque areas. A similar star - count ing 
technique, but with a fixed cell spacing, has been used by lDobashi et al. (Hooi) to 
make an atlas of dark clouds using DSS-I data. These authors identify 2448 dark 
clouds toward the galactic plane (6 < 40°) and present maps with angular reso- 
lutions of 6' and 18' that are publicly available (http : //darkclouds . term. jp). 

Near infrared observations can be used to extend the optical extinction mea- 
surements to the most opaque regions of clouds thanks to the wavelength depen- 
dence of the dust absorption. The recent development of large-format infrared 
arrays has made it possible to apply this technique to increasingly large regions of 
the sky, and thus complement the work carried out at optical wavelengths. In ad- 
dition to applying the star-count method, which provides an estimate of the aver- 
age extinction inside a grid cell, multi-band NIR observations can be used to pro- 
vide a direct determination of the extinction towards individual background stars. 
This determination uses the small range of NIR colors spanned by the background 
stars (mostly giants), which for practical purposes and with little e rror can be 



assumed constant. The Near Infrared Color Excess (NICE) method of lLada et al. 



uses H and K observations of background stars together with an estimate 
of their intrinsic (H-K) color from a nearby unobscured con trol field to derive 
Ay for each observed star. As an application of this method, Lada et al. ( 19941 ) 



produced an accurate map of the e xtinction towards the IC 51 46 dark cloud with 
a resolution of 1.5 arcminute (also Lada. Alves &: Lada 19991 ) . An improvement 
of thi s technique is the NICER (NICE Revisited) method of iLombardi &: Alves 
(|200ll ). which uses observations in additional NIR bands to reduce the variance 
of the Ay esti mate. By applying thi s method to the Orion data of the all-sky 
2MASS survey (ISkrutskie et al.ll2006l ). these authors show that it is possible to 
derive large-scale cloud structure with 5' resolution to a limiting extinction of 
Ay ~ 0.5 (better or comparable to optical star-count work and below the thresh- 
old for molecular CO formation). Future use of this technique with 8m class 
telescopes is expected to provide accurate extinction maps with two orders of 
magn itude dynamical range and a resolution better than 10" (ILombardi &: Alves 

200 ih . 
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2.1.2 RADIO DATA In parallel with the technical advances in NIR cam- 
eras, millimeter and submillimeter-wavelength bolometer arrays have grown in 
size and sensitivity over the last decade. Arrays like SCUBA on the JCMT, 
MAMBO on the IRAM 30m telescope, and Bolocam on the CSO have made it 
possi ble to map systematically the thermal emission of cold dust from dark c louds 
(e.s. iMotte. Andre fc N^ll998l : [fohnstone et aP 1200(11 : lE^ch et al.l [2nnfil V In 
contrast with their NIR counterpart, however, millimeter (mm) and submillime- 
ter (submm) observations of dark clouds are sensitive to a narrower range of dust 
temperatures and suffer from severe instrumental limitations. The need to chop 
the telescope between sky and reference positions filters out part of the large- 
scale emission, which is usually weak enough to be close to the detection limit 
of the array. This means that bolometer observations miss most of the extended 
emission from a cloud while they recover its full small-scale structure. The future 
SCUBA-2 array at the JCMT, with its expected thousand-fold increase in map- 
ping speed and high sensitivity for extended emission, promises to revolutionize 
the field of dust continuum mapping of dark clouds. 

Molecular-line observations, on the other hand, provide information about the 
gas component of clouds, in particular about their velocity structure, density, 
temperature, and chemical composition. These data are complementary to the 
dust extinction and emission measurements, and they are required to under- 
stand, among other issues, the equilibrium status and the chemical evolution of 
dark clouds. Progress in mm-line receiver technology over the last ten years has 
also improved our ability to map large areas of the sky with increasing angular 
and spectral detail. Low resolution (~ 10') CO ir naging of the entire Milky Way 
has b een carried out with the CfA 1.2m telescope (iDame. Hartmann &: ThaddeusI 
200 ih . and systematic mapping of selected clouds at higher resolution (~ 2') has 
been performed in CO isq topologues with the NANTEN/Nag oya 4m telescop e 
fe.g.. lTachihara et aLlliooi ) and the KOSMA 3m telescope (;e.g.. lSun et"aDl200fil ). 
The 32-beam SEQUOIA array (and its predecessor, the 15-beam QUARRY ar- 
ray), operating at 3mm on the FCRAO 14m telescope, have been used to make 
megapixel i mages with sub-arcin inute resolution of nearby clouds like Taurus 
(Fig. 1 andlG oldsmithet_alJl2QQ5|) and furth er out regions like the Galactic Ring 



( Jackson et al.. .2006. ) and the outer Galaxy ( Heyer et al. 19981 ). This array has 



also provided line data for the Coordinated Molecular Probe Line Extinction and 
Thermal Emission (COMPLETE) project (IPidge et al.ll200()l ). which has the goal 
of combining the different IR and radio techniques of dark cloud mapping to study 
a series of dark clouds common to the Spitzer Space Te lescope Legacy Pr ogram 
"Prom Molecular Cores to Planet Forming Disks" (c2d. lEvans et al.ll2003l ). The 
potential of heterodyne arrays, however, has not been fully exploited yet because 
the reduced number of pixels in the receivers limits the large scale observations 
to the bright lines of CO isotopologues. Further progress in this technology is 
still crucial to our understanding of the large-scale physics and chemistry of dark 
clouds. 



2.2 Observational Properties of Dark Clouds 

The dark cloud data obtained using the new instrumentation complem ents a large 
body o f previous work summarized i n pre vious reviews like those bv iGoldsmith 
(|l987l l. iLada. Strom k Mverd (|l993l l. and iMverd (|l995l l. These new data have 
helped greatly to refine our understanding of the dark cloud properties, although. 
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as we will see, a number of important uncertainties still remain. In the following 
sections, we review some of the basic properties of clouds with emphasis on 
recently clarified issues and general characteristics from an observational point of 
view, concluding with a brief review of the issues whose solution remains pending. 
Because of space limitations, we cannot do justice to all the excellent studies of 
individual clouds carried out with the optical, IR, and radio instrumentation 
discussed before. 



2.3 Shapes 



A simple inspection of dark cloud images obtained by any of the methods dis- 
cussed previously reveals that clouds come in a variety of sizes and shapes. In 
general, dark clouds have highly irregular edges, and their overall appearance is 
filamentary and often wind-blown. The pr esence of l ong, well-defined filaments 
was emphasized already a century ago by Barnard (1903) when he noted, de- 
scribing the Taurus plate of Fig. 1, that "among the most surprising things in 
connection with these nebula-filled holes are the vacant lanes that so frequently 
run from them for great distances." Indeed, some of the filaments in Taurus 
can be followed for more than 4 degrees or 10 pc, and similarly thin and long 
structures can be seen in r nany other clouds like Ophiu chus, Lupus, and Orion 



both in optical p l ates (e.g., Schnei der & Elmegreen 19791 ) and radio images (e.g. 



Ballv et al.lll987l : iJohnstone fc Ballv, 19991 ^. 

In many clouds including Taurus and Ophiuchus, the length of some filaments 
is comparable to the full extent of the cloud. Typically, a cloud contains two or 
three long filaments that are either parallel or converge at a low a ngle in a massive 



cond e nsation that often contains an active cluster-forming site (jTachihara et al. 



2OO2I : iBurkert Hartmamil l2004l ) . The veloc ity field of some filaments, in addl 



tion, seems rather coherent (see lLorenlll989l . for a study of Ophiuchus). This 



combination of spatial length and velocity coherence for some of the filaments 
suggests that their presence is intrinsic to the cloud structure and not the result 
of later evolutionary factors like star-formation activity, which would produce a 
more chaotic, small-scale mass distribution. Dark clouds, therefore, seem to be 
born with a filamentary distribution of material that extends over a number of 
parsecs. As clouds evolve and form stars, the pr oducts of star fo rmation inherit 
the filamentary distribution of the parental gas (|Hartmannll2002l ;). 



2.4 Mass Distribution 

Maps of dark clouds provide more than just information on shapes. Measure- 
ments of dust extinction or gas column density at each cloud position can be 
used to estimate the amount of material under different physical conditions, and 
in particu l ar at different densities. In his extinction study of 24 dark clouds, 
Cambresvl (|l999l ;i finds that most of them present a similar power-law relation 



between the mass enclosed in an iso-extinction contour and the extinction, al- 
though the Ay range of this optical study is restricted to Ay ~ 1-5. Despite this 
limitation, the measurements illustrate how most of the material in a cloud lies 
at relatively low extinctions, and it must therefore r eside in the form of low den - 
sity gas. Using more sensitive NIR extinction data, Alves. Lada Lada ( 19991 ) 



find that only 25-30% of the total m ass in the IC 5146 and L977 cloud lies at 
Ay > 10, while Cambresv et al. ( 20021 ) finds that this percentage is slightly more 



Cold Dark Clouds 



7 



than 10% for the North America Nebula. A much smaller value (~ 1%) has been 
recently reported for the Pipe Neb ula, also from IR extinction m easurements 
( Lombardi. Alves &: Ladair2006l . also lLada. Alves Lombardill2007l ). Similar low 
fractions of dense material are indicate d by large-scale maps of dust eni ission 
in clouds ac t ively forming stars. Johnstone. Di Francesco &: Kirlj ( 2004 ) and 
Young et al. find that only between 1 and 2% of the mass in the Ophiuchus 

cloud is associated with dense, submm continuum emission, and that this emis- 
sion is only detected towards regions of extinction larger than 10-15 mag. This is 
illustrated in Figur e 2, which sh o ws th e IR-extinction and submm-emission maps 
of Ophiuchus from lYoung et al.l ( 20061 ). 



A related measure of how the mass of a cloud is distributed amo ng differ- 
ent regimes is the statistics of extinction values. Ridge et al. ( 20061 ) find log- 
normal distributions for the number of pixels as a function of extinction in 
both Perseus and Ophiuchus, a distribution that is i n agreement with the ex- 
pectation from turbulent model s of cloud formation (IVa zquez-S emadenil Il994l : 
Ostriker. Stone &: Gammi3l200ll ) . iLombardi. Alves &: Lada (2006,) , in their study 
of the Pipe Nebula, find that while a lognormal distribution provides an ap- 
proximate general fit to the data, the observed histogram has significantly more 
structure that can only be reproduced with additional components. 



2.5 Velocity Structure 

The presence of supersonic motions in dark clouds has been known since the first 
molecular-line studies revealed non-gaussian line shapes and velocity differences 
across clouds of the order of 1 km s~^ or more (sound speed is 0.2 km s~^ for 10 
K gas). These motions are inferred to be turbul ent from the lack of systematic 
patterns like rotation, expansion, or infall (e.g., iHeiles fc Kat j 19761 ). and from 
the existence of a systematic power-law relation between the velocity dis persion 
inside a cloud a and its physical size L (a(km 8"^)= 1.10 L(pc)°-3^, iLarsonl 
[l981. ). which is reminiscent of the classical Kolmogorov law of (incompressible) 
turbulence. Understanding and characterizing these turbulent motions, which 
affect not just dark clouds but the whole interstellar medium at almost every 
scale, is a formidable task. In contrast with their laboratory counterparts, turbu- 
lent motions in the interstellar medium are supersonic, compressible, and likely 
of magnetohydrodynamical nature. A rich field of both analytic and numerical 
work has been developed to tackle this problem, and a review of its results ex- 
ce eds the scope of this paper . We refer the reader to the excellent recent review 
by lElmegreen &: Sea U (120041 ^ for an in-depth view of interstellar turbulence, es- 
pecially from a theoretical point of view, and for an extensive list of references. 
Here we will concentrate on the main observational results of recent studies of 
dark cloud kinematics. 

When studying the kinematics of cloud gas, it is necessary to take into account 
its dependence on density. Models of supersonic turbulence show that the highest 
density gas moves slowe r than the low den sity material, as it arises from the 
convergence of flows (e.g., Padoan et al.ll200"l ). and this velocity segregation seems 



to be present in clouds. The densest regime in a cloud corresponds to the dense 
cores, which a re known to h ave very low levels of turbulence and subsonic internal 
motions (e.g., Myers 19831 ). Dense cores seem therefore to lie at the bottom of 
the hierarchy of cloud kinemat ics, and to represent s tructures dominated by the 
dissipation of turbulence (e.g. Goodman et al. 19981 ). In § [3] we will present a 
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more complete summary of the physical properties of dense cores, and defer the 
discussion of dense gas kinematics to that section. 

The kinematics of the low density gas in dark clouds a t the small scales typica l 
of dense cores (~ 0.1 pc) has been studied in detail by gar one et al.l 
These authors selected 3 regions in nearby clouds and mapped them in different 
CO isotopologues over several arcminutes (tenths of a parsec for the clouds dis- 
tance) with the IRAM 30m telescope. Each selected region contains a starless 
dense core, although the CO emission is dominated by the more diffuse material 
because of optical depth and depletion effects (see ^ 14.21). Ev en at their limit 
of angular resolution (about 1700 AU), Falgarone et al. ( 19981 ) find unresolved 
velocity structure, especially at the highest speeds. This fastest gas appears 
as prominent wings in non-Gaussian CO spectra with velocities of about 1 km 
s~^ and has an extremely filamentary distribution with aspect ratios larger than 
about 5. In contrast with the dense (10^ cm~^), subsonic, gas in a core, the 
fast CO-emitting gas (with a density of a few 10^ cm~^) presents large shears 
and velocity dispersions at the smallest resolvable spatial scales. The small re- 
gions of high shear may represent locations of enhanced dissipation of turbulence 
(|Petv Falgaronell2003l 'l. 

The large-scale distribution of turbulence in clouds has been studied by a num- 
ber of authors using d i ffereii t techniques. In addition to the classical linewidth- 
size relation of Larson ( 198ll ). methods used to an alyze cloud kinematics includ e 
the statistics of the centroid velocity flu ctuations (Miesch. Scalo &: Bally 19991 ). 



the spectral correlation fun ction (SCF, Rosolowsky et al. 19991 ). the principal 



component an alysis (PCA, Hever &: SchloerbT 1997I ). and the velocity channel 
analysis (VCA, iLazarian &: Pogosyan 200d ). A common goal of these methods is 
to identify and extract from maps of molecular-line emission a number of statis- 
tical properties of the velocity field that can be compared with the expectations 
from analytic theory (like Kolmogorov's law) or from num erical simulations of hy- 
drodynamic and magnetohydrqdynam ic turbulence (e.g.. lOssenkopf &: Mac Low 
2OO2I : iBallesteros-Paredes et al.ll2007f ). Two main issues addressed by these stud- 



ies are the way in which the turbulent energy is distributed between the different 
spatial scales (as described by the energy spectrum E(k)) an d the nature of the 
driving agent (for more details, see lElmegreen fc Scalol[2{)o3 ^. 

A recent example of the above type of work is the principal c omponent anal- 
ysis o f 23 fields of the FCRAO CO Ou ter Galaxy Survey by iBrunt Heyer 
(j2002l ). Using these and additional data, iHeyer Sz BruntI (j2004l ) have found an 
almost universal power-law form for the energy spectrum of clouds, with very 
little cloud-to-cloud scatter (10-20%) not only in the power law index but in the 
normalization coefficient. Such little scatter would help explai n the ernergen ce of 
Larson's law when combining data from an ensemble of clouds ( LarsonI 1981 ). By 
combining numerical simulations of supersonic turbulence with an analysis of the 
velocity line centroid, the linewidth-size r e lation , and the PCA determinations of 
turbulence properties. Brunt Mac Lowl ( 2004 ) suggest that observations of line 
emission from clouds are best fitted by compressible, shock-dominated (Burgers) 
turbulence. They also suggest that either clouds are recently formed or they 
are continually driven on large scales. These suggestive results still need confir- 
mation from the systematic observation of a large number of nearby clouds and 
from the realistic modeling of the combined effect of cloud kinematics, chemistry, 
and radiative transfer. The extreme complexity of turbulent motions in clouds, 
the diverse mechanisms likely involved in their driving and dissipation, and the 



Cold Dark Clouds 



9 



challenge of obtaining high-quality molecular-line data of whole clouds will surely 
keep challenging theorists and observers alike for a number of years. 



2.6 Internal Structure 

A common feature of dark clouds is that they show a partly hierarchical structure, 
with smaller subunits appearing within large ones when observed with increas- 
ing spatial resolution. To characterize this structure, two different approaches 
have been generally followed, depending on whether the discrete or the contin- 
uous nature of the structure is emphasized. In the first approach, the cloud is 
assumed to be composed of subunits, referred as clumps and defined as coher- 
ent regions in position- velocity space that may contain significant substructure 
( Williams. Blitz &: McKej 2000l ) (see below and Table 1 for typical clump pa- 
rameters) . These clumps are identified from the data using an automatic clump- 
finding algorithm that simplifies the process and avoids hum an bias. The most 
popul ar algorithms in use are G AUSSCLUMPS, developed b vlStutzki k Giistenl 
(|l990l l. and CLUMPFIND from IWilhams. de Geus Blit d ^mi). These algo- 
rithms follow different approaches to identify and characterize clumps, but they 
generally agree in the derived properties, es pecially in the intermediate and high 
mass end ( Williams, de Geus fc Blit j|l994l ). When they are apphed, usually to 
CO isotopologue data, the spectrum of clump properties is found to follow a power 
law in the range for which the data are complete. The mass of the clumps, for 
example, commonly presents a distribution of the form dN/dM ~ M~", w here a 
lies in a narrow range betw een 1.4 and 1.8 for both dark clouds and GMCs (jBlitd 
19931 : iKramer et al.| [T998l ). As a reference, and in this form, we note that the 
i nitial mass fun ction of stars follows a steeper power law with a slope of a = 2.35 
(jSalpeteiill955l ). The clump distribution, therefore, contains most of the mass in 
massive clumps, while the stellar distribution has most of the mass in low-mass 
objects. 

The second approach to characterize cloud structure assumes that the cloud 
is self-similar, at least over some range of scales, and applies the concepts of 
fractal geometry. In this geometry, the most characteristic parameter of an ob- 
ject is it s fractal (or Haus dorff) dimension, which is in general a non integer 
number ( MandelbrotI Il983 ) . The boundary of a cloud, for example, has a frac- 



tal dimension Dp that can be estimated from the relation between the enclosed 



area {A) and the perimeter length (P): P ~ ( Love joy 19821 ). For sim- 
ple curves in a plane, like circles or ellipses, this area-perimeter relation gives 
the expected value of Dp = 1 for the dimension of a line. For highly frag- 
mented or convoluted curves, the expression gives a value that is large r than 
1 and that approaches 2 as the curve fills a larger fraction of the plane. iBeech 
(|l987l ) used the above area-perimeter relation for the optical boundary of 24 dark 
clouds and derived an average boundary dimension of 1.4, which suggests that 
clouds are fractal. Using a similar n i ethod for the iso-contours of a series of CO 
maps, iFakarone. Phillips fc Walked (Il99lh have derived a bounday dimension 



Dp = 1.36 lb 0.02 for the Taurus molecular cloud, almost the same value mea- 



sured by Scalo ( 199d ) from IRAS maps. Similar values for the f ractal dimension 
have been obtained for the more diffuse cirrus clouds {Dp = 1.3, Bazell fc Desertl 
19881 ) and even for the atmospheric clouds on Earth {Dp = 1.35. lLoveiovlll982l ). 
As these dimension estimates use plane-of-the-sky views or maps of the objects, 
they only provide the fractal dimension of the clouds projection. To infer the 
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intrinsic (volu me) fractal dimension, the area-perimeter value is generally in- 
creased by 1 ( Beech 19921 ). although this procedure may only provide a lower 
limit to the true dimension ( Sanchez. Alfaro &: Pere3 20051 ). With this assump- 
tion, the area-perimeter determinations imply a typical cloud fractal dimension 
of around D = 2.4. 

The very different approaches of the clump-finding and fractal analysis would 
seem to suggest that the two methods provide contradictory views of the cloud 
internal structure. This is however not the the power-law behavior found 

in the spectra of clump proper ties already suggests a scale- free s tructure con- 
sisten t with a fr actal geometry ( Elmegreen &: Falgarond 1996 , also Stutzki et al 



19981 ). Indeed, Elmegreen Sz Falgarong (|l996l ) show that it is possible to de- 



rive clump mass distribution functions as those found by the GAUSSCLUMPS 
and CLUMPFIND algorithms from fractal cloud models of (volume) dimension 
D = 2.3 ± 0.3. In this view, the clump-finding algorithms identify the peaks of 
the fractal intensity distribution of a cloud. Due to the threshold imposed in the 
density or column density by the observations, a power-law spectrum of cloud 
sizes results, and this produces a power-law distr ibution of clump masses even if 
the gas density follows a lognormal distribution ( ElmegreenI 2002 ). 

The self-similar behavior revealed by both the power-law clump decomposi- 
tion and the fractal description seems to apply inside a given range of spa- 
tial scales. This range most likely represents the regime dominated by turbu- 
lent r notions, as turbule r ice p roduces self-similar structures in a natural way 
(e.g., Elmegreen &: Scalol 20041 ) . At large scales, the presence of well-defined 
features like filaments or shells indicates a deviation from self similarity that 
is likely caused by the mech anisms responsible for the cloud formation (e.g., 
Elmegreen fc Falgaronjll99fil ). 



At small scales, in the gravity-dominated regime, the self-similar picture is also 
expected to break down. No break point from self similarity has yet been found 
in the lower density gas, and iFalgarone et al. I (Il998l ) report unresolved turbulent 
structures at scales as small as 0.02 pc from their study of the highest velocity 
gas in the Taurus molecular cloud ( also IPetv &: Falgarone I2OO3I ). For the dense 
star-forming gas, on the other hand, iLarsonl ^iml ) has inferred a break down of 
the self-similar behavior in the form of a discontinuity in the slope of the surface 
density of companions for young stars in Taurus. This break point divides the 
binary regime from the clustering regime, and it occurs approximately at a radius 
of 0.04 pc, which coincides with the Jeans length of the cloud. Such coincidence 
suggests that for scales smaller than this length, thermal pressure provides the 
dominant support against gravity^ whil e for larger scales, tu rbulent and ni agnetic 
pressures are important ( Larson! 19951 ). Followup work by Simon (|l997l ) shows 
similar break points in the density of stellar companions for Ophiuchus and the 
Orion Trapezium, although the scale-size of the break point does not follow the 
expected Jeans mass. A larger value f or the break point in th e self-similar be- 
havior of Taurus has been proposed by Blitz &: Williams ( 199?! ) based on ^^CO 
observations, although optical depth effects and molecular freeze out (section 
§ 14. 2p may have affected this estimate. 

There is further evidence for a change in the structure of dark clouds at small 
scales. As mentioned before, the clump mass spectrum found from large-scale CO 
observations is relatively flat, and follows a power law of the form dN/dM ~ 
with a = 1.4-1.8. The mass spectrum of millimeter/submillimeter dust con- 
tinuum peaks, on the other hand, presents a steeper slope of a =2.0-2.5 for 



Cold Dark Clouds 



11 



Table 1: Properties of Dark Clouds, Clumps, and Cores 



Mass (Mq) 
Size (pc) 

Mean density (cm~^) 
Velocity extent (km s^ 
Crossing time (Myr) 
Gas temperature (K) 
Examples 



Clouds" 
10^- 



-10^ 
2-15 
50-500 
2-5 
2-4 
10 

Taurus, Oph, 
Musca 



Clumps 



Cores'^ 



50-500 
0.3-3 
10^-10^ 
0.3-3 
« 1 
10-20 
B213, L1709 



0.5-5 
0.03-0.2 
10^-10^ 
0.1-0.3 
0.5-1 
8-12 
L1544, L1498, 
B68 



° Cloud masses and sizes from the extinction maps bv iCambres^ (|l999l ). 

velocities and temperatures from ind ividual clo ud CO studies 

^ Clump properties from Loren ( 19891 ) (^^CO data) and 



Wilhams. de Geus Blitzj (119941) fCO data) 

Core properties fro m Jiiina. Mvers &: Adamsl ( 19991 ). ICaselli et al.l ( 2002al ). 
Motte. Andre Sz Neri ( 19981 ). and individual studies using NH3 and N2H+ 



mass e s larger than about 1 ( Motte. Andre &: Neril Il998l : iTesti Sargent! 
19981 : Johnstone et al. 2000l . 200ll ). This spectrum, whic h at least for the case 
of Ophiuchus flattens to ab out a = 1.5 for lower masses ( Motte. Andre fc Neril 



19981 : 1 Johnstone et al.ll2000l). mimics the shap e of the initial mass function (IMF) 
of stars (jSalpete il ll95,5i r iMiller Scald Il979l ). in contrast with the clump mass 
spectrum. Such similarity to the IMF suggests that stellar masses may be de- 
termined by the same process that fragments the molecular gas at the smallest 
scales. 

The above discontinuity in the slope a still needs confirmation with a single 
observational technique, as each side of the break point has been observed with a 
different tracer (CO for large scales and mm/submm continuum for small scales). 
If confirmed, the discontinuity will suggest a change in the physics of the cloud 
gas at the tenth of a parsec scale, and this can provide the basis for an empirical 



descri ption of the different levels of cloud structure. IWilliams. Blitz McKee 
( 2OO0I ) have recently proposed one such a categorization of structure in terms 
of clouds, clumps, and cores, with clumps being defined (as before) by their 
velocity coherence and cores being defined as gravitationally bound, single peaked 
regions out of which individual stars (or sim ple stellar systenis) forni. Alte rnative 
definitions have already been proposed (e.g.. lGoldsmithlll987l : lMverslll995l ) . and it 
is also possible that any description of cloud structure using o nly a few element s 
is too simple to capture the continuous properties of cloud gas (jRodrigue 3 l200,5l ). 
Still, the need to name and characterize the different levels in the hierarchical 
structure of cloud s makes it necessary to have a wel l accepted terminology, and the 
categorization by Williams. Blitz &: McKee ( 2000l ) matches the spirit of current 
usage. Thus, using this convention, we present in Table [Da summary of the main 
physical parameters of clouds, clumps, and cores. We stress that these units 
are still loosely defined and that their properties may be sensitive to the tracer 
used in the measurement and can vary with cloud. Future improvements in our 
understanding of cloud chemistry and structure will help further refine the above 
description. 
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2.7 Magnetic Field 



The magnetic field of a cloud is probably its most difficult property to measure. 
The line of sight strength of the field can only be directly determined observing 
the Zeeman splitting of line transitions, and its plane of the sky direction can be 
estimated via polarization measurements of background stars, dust emission, or 
spectral lines. All these measurements require observations that combine stability 
and high signal t o noise, and are therefore difficult and time consuming (see 
Heiles et al.lll998l for an in depth review of observational techniques of magnetic 
field measurements). Although the magnetic field only acts directly on charged 
particles (electrons, ions, and charged dust grains), its presence can be felt by 
the neutral material through collisions. Under most dark cloud conditions, this 
ion-neutral coupling is highly efficient and , except for the de nsest regions, the 
field is expected to be frozen to the gas (see iMcKee et alll99.'j for a review of the 
basic magnetic field theory). Because of this, the ability of the magnetic field to 
counteract the action of self-gravity can be critical to the equilibrium balance of 
dark clouds. 

The most straightforward way to map the large-scale orientation of the mag- 
netic field in a cloud is to measure the polarization of light from background 
stars, which results from the dic hroic extinction by aligned, elongated dust grains 
(e.g., Davis &: Greenstein 195ll ). Large-scale m a ps of polar ization for t h e Tau - 
rus cloud have been p resented by iMoneti et al.l ( 19841 ) and iHever et al.l ( 19871 ). 
iGoodman et al. (1990) have produced maps of the optical polarization for Perseus, 
Taurus, and Ophiuchus using a combination of their own measurements with pre- 
vious data. From these maps, they find that in Taurus and Ophiuchus the pat- 
tern of polarization is highly regular over scales of about 10 pc, while the Perseus 
measurements seem to suffer from the superposition of two components along the 
line of sight. In Taurus and Ophiuchus, the large scale filaments are neither com- 
pletely parallel nor completely perpendicular to the global polarization pattern, 
sugges ting that the magnetic field does not dominate the cloud structure on large 
scales ( Goodman et al. 19901 ). Optical polarization data from other clouds offer a 
mixed picture. Clouds gen erally present fields with a well-defined mean direction 
but significant dispersion (jMvers &: GoodmanI Il99ll ). and cases where the field 
is parallel or perpendicular to a filament exist (e.g., the long filamentary Musca 
cloud presents an ordered fi eld perpendicular to the long axis of the filament, see 
Perevra &: Magalhaeslliooi ). 



While optical polarization measurements are necessarily limited to the low- 
extinction parts of a cloud, so the background stars are still bright enough to have 
their polarization measured, IR polarization data could in principle sample more 
opaque regions and therefore provide a deeper view of the inner magnetic field. 
Unfortunately, this seems not to be always the case. For the Taurus cloud, for 
example, there is evidence that the dust loses its polarizing power at depths larger 
than about Ay = 1.3, so IR pola rization measurements offer little impro vement 
over optical data in this region ( Goodman et al. 19951 : Arce et aD 19981 ). This 
loss of polarizing power probably results from a combination of a change in the 
alignment properties of grains and a change in their optical properties because 
of coagulati on and mantle growth, although the exact mechanism is not yet well 
understood ( Lazarian. Goodman k Myers[ll997l : Twhittet et al.ll200ll ). 

Better sampling of the innermost cloud regions is expected from polariza- 
tion measurements of the millimeter /submillimeter emission from dust, although 
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they also suffer often from drops in polarizing power towards the densest ar- 
eas ( Matthews. Wilson &: Fieg3 200 ll ). Because of limitations due to sensitivity 
and spatial filtering, large scale studies of the dust polarization have concen- 
trated on the brightest objects like the Orion filament ( Matthews &: Wilson! 2000 : 



Houde et alJ 12004) instead of the more nearby and colder dark clouds (see 
for submm-polarization studies of dense cores). As in the case of the optical/IR 
polarization measurements, the limited submm measurements do not show a very 
strong correlation between the magnetic field direction and the mass distribution 
in the clouds (e.g., iHoude et al. 2004). 

Zeeman effect observations of dark clouds provide an estimate of the strength 
of the (line of sight component of the) magnetic field. The high signal-to-noise 
required by these measurements and the low spatial resolution achieved, because 
of the low frequency of the transitions, means that no large-scale maps of Zee- 
man effect measurements are currently available. Observations are commonly 
restricted to single pointings, which offer only a limited sampling of the magnetic 
field strength in clouds. In addition, and due to the need for strong emission, 
the pointings tend to be selected towards the densest regions. Still, coming up 
with meaningf ul detections is challe nging, and the OH Zeeman observations of 
dark clouds bv lCrutcher et al" ( 19931 ) produced one clear detection out of 12 po- 
sitions observed (n iostly in Taurus and Ophiuchus). These data, together with 
the observations of iTroland et al. ( 19961 ) suggest typical magnetic field strengths 
of order of 10 /iG or less. 

A main issue for the studies of the magnetic field strength in cloud gas is the 
importance of the field in the dynamics and balance of the cloud. A convenient 
parameter to characterize this importance is the observed ratio between the mass 
and the magnetic flux normalized to the critical value aG~ ^^'^, where G is the 
gravitational constant and a ~ 0.13 fsee lMcKee et al.lll993l for a full discussion 
and equations). If the observed mass-to- flux ratio exceeds the critical value, the 
cloud is said to be supercritical and the magnetic field cannot prevent collapse. 
If the ratio is lower than the critical value, the cloud is called subcritical and 
gravity is stabilized by the magnetic field. A detailed compilation and analysis 
of the currently av ailable data on m agnetic field strength in molecular clouds has 
been pres ented bv ICrutched (I1999I') . and an update of this work appears in the 
review bv iHeiles fc Cr'IItS I2OO5I ) . According to this analysis, the magnetic 
field strength is just at the level of being critical, i.e., at the boundary between 
being dynamically important and not. Even in this best-to-date analysis, there 
are important caveats due to the large number of non-detections and the need 
to correct geometrically the mass-to- fiux ratio (inultipl ying by 2 in case of a 
spher e or by 3 in case of a sheet, see iBourke et al.l I2OOII and iHeiles fc Crutcher 
2OO5I I. The values, in addition, are likely more representative of the densest 



regions than of the extended cloud gas although chemical effects may further 
complicate the interpretation (^ 14. 2p . Thus, despite the enormous observational 
effort, the dynamical importance of the magnetic field in clouds, especially at the 
large scales, remains elusive. The current best guess is that the magnetic field 
contributes in a non-negligible way to the energy balance of molecular clouds. 
Whether it is a dominant player or a second order effect cannot be decided yet, 
and because of this, our understanding of the global equilibrium of dark clouds 
is still incomplete. 
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2.8 Equilibrium State and Star Formation 



The equilibrium state of dark clouds and the way they collapse under gravity to 
form stars is probably the most controversial issue related to their nature. The 
uncertainties in the magnetic field strength and the importance of the turbu- 
lent mo tions, together with a possible revision of the l ifetimes of the molecular 
gas re.g. lHartmann. Ballesteros-Paredes fc BerginlB ) , have resulted in two op- 
posed views of the global state of clouds. Briefly, one view holds that clouds are 
close to equilibrium and that their evolution toward star formation is approxi- 
mately quasistatic. The other view defends that clouds are dynamic objects that 
evolve and form stars in a crossing time. The limited space of this review can- 
not make justice to the number of issues involved in this controversy or to the 
detailed position of each camp, which consists of a number of authors working 
often independently. Here we will simply review the main contentious points and 
mention a number of relevant references as an introduction to the topic. We refer 
the reader to the original papers and the reviews mentioned below for an in-depth 
view of the controversy. 

The quasistatic view holds that clouds are objects close to equilibrium, due 
to their relatively long lives (at least 10 Myr according to the recent estimate 
by Mouschovias. Tassis &: Kunz 2006 . see also Blitz &: Shu 198d) and their beiri g 
gravitationally bound and close to virialized (e.g., Larsoni " Il98ll : lMcKeel[T999l l. 
In this view, the equilibrium against self gravity is provided by the magnetic 
field, which has a twofold contribution toward stability. If the static compo- 
nent of the magnetic field is strong enough to make the cloud subcritical, grav- 
itational forces cannot overcome magnetic fo r ces as long as the field re mains 
frozen into the matter ( Mestel &: Spitzer 19561 : Nakano &: Nakamura 19781 ). Hy- 
dromagnetic waves, in addition, can provide additional su pport and contribute 
to the supersonic inotion s observed in molecular clouds (jArons &: Max! I1975I : 
Gammie &: Ostriker 19961 ). In this magnetically-dominated picture, molecular 
clouds can only evolve under gravity and form stars through the process of "am- 
bipolar diffusion," by which neutrals drift past the ions and the magnetic field, 
which remain frozen to each other. Through the action of ambipolar diffusion, 
cores of dense gas form by gravitational contraction out of the initially subcritical 
medium. When these cores have a ccumulated enough mass, they become super- 



J f II 1 

critical and collapse to for m stars (jMestel &: Spitzeiill956l : IShu. Adams &: Lizand 



19871 : iMouschovias fc Ciole k 199^). Under typical cloud conditions, ambipolar 
diffusion is slow (several cloud dynamical times, see § I4.3.2p . and this slowness 
makes cloud evolution prior to star formation occur quasistaticly. It also makes 
star formation a rather inefficient p rocess, in agreement with the low rate of 
stellar birth observed in the Galaxy (IZuckerman Evanslll974l ). 

The more recent and opposite view of cloud evolution and star formation em- 
phasizes the role of supersonic turbulence and lack of equilibrium. Numerical sim- 
ulations have shown that magnetic turbulence decays in about a dynamical time 
(iMac Low et al.lll998l : IStone. Ostriker fc Gami^ll998l : IPadoan fc Nordlundlll999l ) , 
so hydromagnetic waves cannot provide support against gravity without a source 
of continuous replenishment. If in addition fields are weak enough so molecular 
clouds are supercritical (e.g., Nakanol 19981 ). cloud evo lution and star f ormation 
becomes a fast process that occurs in a crossing time ( ElmegreenI 2000l ). In this 
picture, clouds form from convergent flows, evolve, and dissipate rapidly, without 
ever reaching a state of equilibrium, and with a typical molecular cloud lifetime of 
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3-5 M yr ( Ballesteros-Paredes. Vazguez-Semadeni &: Scald 19991 : Hartmann^ Ballesteros-Paredes Bergin 



200ll : Vazguez-Semadeni. Ballesteros-Paredes fc Klessen 20031 : Hartmann 20031 ) . 
Molecular clouds are defined observationally by the detection of CO emission (as 
opposed to H2 emission) and this lifetime refers to the age of the CO emitting 
cloud, which is associated with star formation. Chemical models su ggest that H2 
formation is a pre-requisite to CO formation ( Bergin et al. 20041 ) . and lifetime 
estimates do not include any earlier phase when gas is predominantly molecular 
(e.g. H2) but CO has not formed in sufficient quantity for detectable emission. 
Factors such as ram pressure, molecule and dust shielding, a.nd gas-grain physics 
influence the timescales for both H2 and CO formation (jKoyama fc Inutsuka 



2OOOI : [ Bergin et al.l [20041 ) . But in any case, it is clear that substantial evolu- 
tion (~ 10 — 20 Myr) occ urs in a pre-CO p hase where the cold H I and H2 



would be difficult to detect ( Bergin et al. 20041 ) . Formation of cloud structure in 



this scenario is regulated by the stochastic action of turbulence, which produces 
strong density perturbations through its fluctuating velocity fiel d. Regions of 



dense gas form at the stagnation point of two convergent flows (iPadoan et al 
2OOII I. although such a transient structure often disperses when the flows fade. 



If the initial compression is strong enough, it may decouple from the medium 
and produce a dense core that subsequently collapses to form stars. This process 
of creating structure through the interplay betwe en gravity and turbulence is of- 



Klessen &: Ballesteros-Paredes 



ten r eferred as gravoturbulent fragmentation (e.g 

2OO4I ) an d lies at the heart of the dynami c picture of s tar formation (see the re- 
views bv iMac Low fc KlessenI l2004l . iLarson .2003 . and iBallesteros-Paredes et al 
20071 ) 



The controversy between the fast and slow modes of cloud evolution and star 
formation is still ongoing and far from settled. Attempts to determine cloud life- 
times using stellar ages, for exam ple, have not only not clarified the issue but have 
led to their own controversy (e.g., Palla &: Stahlerl 20021 : Hartmann 20031 ). Chem- 
ical analysis of the gas can also provide an estimate of cloud timescales. Obser- 
vations of cloud cores, for example, sugge st compositions that are far from equilib- 
riurn and ages that are ^ 3x10^ yrs (e.g. Pratap et al.lll997 ; van Dishoeck and Blake 
19981 ). although this chemical age might not reflect the true cloud age but rather 
refer to the last time the chemistry was reset by some dynamical event. Studies 
of c old H I in dark clouds, on the other hand, suggest longer timescales > 3 x 10^ 
yrs ( Goldsmith &: Li 20051 ). Measurements of the magnetic field, as mentioned 
before, remain ambiguous beca use the da.ta po pulate the boundary between the 
critical and subcritical regimes (|Crutchei]ll999l l and this is often interpreted as 
proof of the magnetic field being dominant or a minor player. 

The above summary of the current controversy illustrates how our understand- 
ing of dark cloud physics is still limited and far from complete. The observations 
and theoretical work carried out over the last decade have clearly revitalized the 
field of cloud studies and challenged many ideas previously assumed correct. Un- 
fortunately, they have not achieved yet a consensus view of how clouds form, 
evolve to form stars, and finally disperse. It is therefore a challenge left for 
the next generation of cloud studies to explain the diversity of observations and 
to develop a global view of clouds, from the large, tens-of-parsecs scale to the 
tenths-of-parsec size of the dense cores, to which we now will turn our attention. 
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3 PHYSICS OF PRE-STELLAR MOLECULAR CORES 



Dense cores are localized density enhancements of the cloud material that have 
been recognized for more than 20 years as the likely sites of low-mass star forma- 
tion in nearby dark clouds like Taurus and Perseus. Initial core studies, based 
on the optical inspection of the Palomar plates together with molecular (espe- 
cially NH3) observations, showed that dense cores have typical sizes of 0.1 pc 
and contain a few solar masses of subson ic material at teinperatures around 10 K 
and average densities of a few 10 cm" -3 (iMvers fc Bensonlll983l : iBenson &: MyergI 
19891 1. The correlation of dense core positions with the location of highly em- 
bedded young stellar objects (YSOs) detected with the IRAS Far-Infrared (FIR) 
satellite soon provided the proof th at some dense cores are currently forming 
stars or have done so very recently ( Beichman et al. 19861 ). The presence of an 



embedded YSO has also led to the classification of cores in two groups: starless 
and star-containing. Strictly speaking, the distinction between the groups is only 
based on the threshold for detection of an embedded YSO with current instru- 
mentati on, which was unti l recently of the order of 0.1 Lq for nearby clouds like 



Taurus ( Mvers et al. 198?! ). The increase in the sensitivity of IR observations 



brought by the Spitzer Space Telescope has led to the discovery of a number of 
so-called very low luminosity objects (VeLLOs) in a few previously- thought star 



less c ores, forcing a revision of the core classification (e.g., L1014, see lYoung et al. 

Still, the distinction between starless and star-containing cores seems a 
fundamental one, as it likely represents the separation between the pre-stellar 
and post-stellar phases of core contraction. The study of starless cores, therefore, 
offers the best opportunity to determine the initial conditions of low-mass stellar 
birth. 

Dark globules (often called Bok globules or Barnard objects) are classically 
defined by their op tical appearance as small, roundish, and dark nebulosities 
(iBok fc R,eillvlll947l l. They constitute a less homogeneous class than the dense 
cores of dark clouds, as they span sizes up to abou t 1 pc an d masses of up 
to almost 10^ Mq, and have lower average densities ( Bok 197?! ). A number of 



small globules, however, shares many characteristics with the dense cores, and 



are a l so the formation sites of individual low-mass stars ( Clemens &: Barvainid 
19881 : lYun fc Clem"m3 Il990l l. These dense globules may have a similar forma- 
tion mechanism as the dense cores, with the difference that globules do not lie 
embedded in a molecular cloud because they have been exposed by an external 
event (such as ionization from a nearby O star, Reipurthl [l983). Like some dense 
cores, some globules seem starless and centrally concentrated, and probably rep- 
resent a phase prior to gravitational collapse. As we will see below, their isolated 
nature provides an additional advantage when studying the internal structure of 
pre-stellar objects ( Alves. Lada &: Lada 200ll ). 

The study of dense cores and globules prior to star formation has progressed 
enormously during the last decade. Earlier work lacked the angular resolution 
needed to study their internal structure, and it therefore concentrated on globa l 
properties like total masse s and average densities (e.g., iBenson &: MversI Il989l : 
Clemens fc: Barvainid[l988l ). The progressive increase in resolution and sensitivity 
of radio and IR observations has finally allowed a detailed view of their interior. 
Through a combination of techniques that trace the gas and dust components, 
a new understanding of their internal structure has started to emerge. Crucial 
to this understanding has been the realization that dense cores and globules. 
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despite their apparent simplicity, have a rather complex chemical composition 
that is responsible for previous contradictions between observations made with 
different tracers. In this section we review some of the recent results concerning 
the physical structure of starless cores and globules, and the following section does 
the same for their chemical properties. For the sake of economy, we will refer 
to both cores and globules simply as cores, and we will only make an explicit 
distinction between the two when necessary. We will also use the term "starless" 
in a broad sense that includes what in the literature is sometimes referred to as 
"pre-stellar." 



3.1 Density 

The density structure of pre-stellar cores is typically estimated through the anal- 
ysis of dust emission or absorption. Three main techniques have been used so 
far: (1) near-infrared extinction mapping of background starlight, (2) mapping 
of mm/submm dust continuum emission, and (3) mapping of dust absorption 
against the bright galactic mid-infrared background emission. Each technique 
uses dust observations to derive a parameter that is proportional to the gas 
column density, and then fits this parameter with a model density profile that 
assumes some simple geometry (usually spherical) . Fig. [3] shows a summary of 
how each technique works, including maps of representative cores (B68, L1544, 
and p Oph D), the main equations used in the analysis, and the radial profiles of 
the dust-derived parameter to be fit with the density model. 

Because of their different assumptions, the above techniques are sensitive in 
different ways to the properties of the dust and to the variation of these properties 
inside the core. Uncertainties in the dust parameters propagate in different man- 
ner and affect differently the final density determination of each technique. The 
NIR extinction method, for example, is independent of the dust temperature, but 
depends on the reddening law (using r^'^ to convert Ak into Ay) and on the 
correlati on between Ay and the gas c olumn density, bot h of which may be vari- 
able (see lRieke fc LebofskylllQSsI and iBohhn et al.l 19781 for typical ISM values). 



The mm/submm dust continuum emission method, on the other hand, depends 
almost linearly on the dust temperature and on the value of the dust emissiv- 
ity (kj^), two parameters that are known to vary systematically as a function of 
core depth (see 13. 2. T] for more details). The mid-infrared absorption method is 
also subject to uncertainties. It is sensitive to fluctuations in the emission from 
polycyclic aromatic hydrocarbons (PAHs), which provides the bright background 
against which the core is observed. It also depends on the somewhat uncertain 
relation between the absorb ing opacity and the h ydrogen column density (see 
Weingartner Sz Draind 2001 and Ragan et al ] l200fii for further details). Given the 



restricted use of each of these methods (due to their difficulty and data require- 
ments), no systematic inter-comparison between them has been carried out so 
far. As shown below, the methods agree in the qualitative characteristics of the 
density profiles, but a quantitative comparison of the inferred densities is still se- 
riously needed. Without this, we can only guess the relative uncertainties of our 
density determination; current wisdom suggests that they are in the factor-of-2 
range. 

A main characteristic of the density profiles derived with the above techniques is 



that t hey require a central flattening. Pioneering work by I Ward-Thompson et al 



(|l994l ) showed that the radial profiles of submillimeter continuum emission for a 
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number of cores could not be fitted using single power-law density distributions 
(assuming a constant dust temperature). At least two power laws were needed to 
fit the data, and the inner power law was always flatter than the outer one, which 
was usually close to r~^. Follow up work with higher sensitivity has confirmed 
these early results for an increasing number of objects, and it has shown that 
as a rule, the density gradient of a core is flatter than r^^ within radii smaller 
than 2500-5000 AU, and that the typical central density of a core is 10^-10^ cm 
(|Andre. Ward-Thompson fc Mottei[l996l : IWard-Thompson. Motte k Andrelll999l ') . 
Although a dust temperature dr op toward the core center is expected from ther- 
mal equilibrium considerations ( Evans et al. 200ll ). the effect seems not strong 
enough to explain alone the flattening found in the continuum data. Observations 
using extinction techniques, which are insensitive to temperature, also find that a 



centra l flattening of the density profile is needed to fit the data. iBacmann et al 



( 2000l ). for example, have modeled MIR (7/im) absorption profiles from a sam- 
ple of nearby starless cores observed with ISOCAM, and they have found that 
a decrease in the slope of the density distribution i s requi red to fit the emission 
in the central 4000-8000 AU. lAlves. Lada Ladal (|200ll l have mapped in de- 
tail the NIR extinction from the B68 globule and achieved an almost perfect fit 
with an equilibrium isothermal (Bonnor-Ebert) density profile that has a central 
flattening within ~5000 AU. Density profiles comparable to those inferred from 
observations are shown in Fig. [H 

A number of density laws with central flattening have be en used to fit the obser- 



vations of dense cores. The dual power law used initially by lWard-Thompson et al 



is illustrative of the need for inner flattening, but presents an unphysical 
discontinuity in the derivative that results ii i a pressure jump. A simple al- 
ternative is the softened power laws used by iTafalla et all (l2002l ). and a more 
phys i cal model is th e truncated isothermal (Bonnor-Ebert, BE) sphere (jEbert 
I955I : Bonnor 19561 ). that often (but not always) provides a good fit to the 
data. Examples of good BE fits include those to the extinction profiles of a 
number of Bok globules (|A ves. Lada &: Ladal I2OOII : iKandori et aP l200,5l ^ and 
dense cores Bacmann et al. Other physical models, on the other hand. 



seem to be ruled out by the se observations. The logotropic density profile of 
iMcLaug hlin &: Pudritj ^dM) is too flat at large ra dii (r ^) to fit most cores 
j^acman n et al.ll2000l ). while the Plummer-like profile ( Whitworth &: Ward-Thompson 
I2OOII ) seems asymptotically too steep (r~^). Less clear is the situation of magnetic 
field-dominated models. These models naturally produce equilibrium configura- 
tions with flattened d ensity profiles and approximately power-law behavior 
at large radius (e. 
Lizano Siuilll98 



MouschoviasI Il976l: iTomisaka. Ikeuchi Sz Nakamural Il988l : 
Bacmann et al. ( 2000l ) find that some of these configura- 



tions provide adequate fits to their data, but that they require magnetic fields 
of the order of 50-150 //G, significantly larger than commonly observed (^ 13. 4p . 
They also predict o blate geometries, which seems to contradict some observations 
(|Mvers et al.lll99lh . 

The high quality of some BE fits raises the question of whether cores are truly 
isothermal spheres solely supported by thermal pressure and confined by an outer 
medium. Several lines of evidence suggest that although thermal pressure is an 
important ingredient in core structure, it is not the only one, and that cores are 
more complex entities than simple BE spheres. In the first place, cores are seldom 
spheri cal, but they appea r in the sky as elliptical objects with axial ratios typically 
of 2:1 (|Mvers et al.lll99lh . Whether cores are prolate, oblate, or triaxial is still a 
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matt er of some debate ( Myers et al. 1991 : RvdenI 1996 : Jones. Basu &: Dubinski 
200ll ) , but the lack of spherical symmetry already suggests either the presence of a 
non symmetrical force component or the lack of perfect equilibrium. In addition, 
the observed density contrast between the core center and its edge often exceeds 
the maximum value of ~ 14 allowed for stable BE spheres (jBacmann et al.l 12000 : 



Evans et al.ll200ll ). and it seems unlikely that such an unstable configuration can 
be realized in nature if thermal pressure is the only balancing force. Finally, the 
pressure required by the BE fits usually exceed the thermal (plus subsonic non- 
thermal) pressure available from the gas by a factor of a few, again i ndicati ng the 
presence of an additional force component or the lack of equilibrium ( Tafalla et al.l 
20041 : iKirk. Ward-Thompson Andrei [2OO5I I. 

The BE- like profiles found when fitting observations must th erefore r esult from 
a more complex phys ics than just thermal hydrostatics. Indeed, MversI ( 20051 ) and 
Kandori et al.l (|2005l ) have shown that these profiles can be reproduced with a va- 
riety of geometries and even under early conditions of gravitational collapse. Even 

models of supersonic turbulence can pro duce density condensations that have in 

proje ction the appearance of BE profiles (jBallesteros-Paredes. Klessen &: Vazquez-Semadeni 
2OO3I ) . but their expected velocity structure is inc onsistent with the prevalence 
of subsonic cores in dark clouds like Taurus (e.g., Jiiina. Myers &: Adamd 1999 . 
Tafalla et al. 2007, in preparation). The ultimate origin of the observed core 
density profiles remains therefore unexplained. 



3.2 Temperature 

The kinetic temperature of the dust and the gas components in a core is reg- 
ulated by the equilibrium between heating and cooling. At typical core den- 
sities (^J^^_cinr£}_;_Jhegasan^^ via collisions 
(e.g. Goldsmith &: Lange 3 ll978l : iBurke fc Hollenbachl[l98^ ) , and the two temper- 
atures are expected to converge; whether they equalize or not depends in part on 
the exact density of the core and on the strength of the external radiation field 
(iGalh. Walmslev fc GoncalvesI I2OO2I ) . ]n this and the next sections we discuss 
recent progress in the determination of the dust and gas temperatures of cores. 
As the two temperatures are determined using different observational techniques 
and are therefore subject to different uncertainties, we separate their discussion. 
We start with the dust. 

3.2.1 Dust Temperature The dust temperature is physically determined 

by the balance between the heating by the interstellar radi ation field (ISRF) and 

the c ooling by thermal radiation from the grains in the FIR (jMathis. Mezger &: Panagia 
19831 ) . To estimate this temperature observationally, color measurements of the 
dust in the submm or FIR are commonly used, sometimes complemented with 
NIR extinction measurements. A major uncertainty in this work is our limited 
understanding of the optical properties of the dust and, in particular, the value 
of the dust emissivity and how it changes as a function of density. Theoretical 
models of dust evolution predict that as the density increases, the dust grains 
become coated with ice ar id coagulate to form fluffy a ggregates, changing their 
emissivity in the process ( Ossenkopf Hennind 1994). Determining dust tem- 
peratures from observations therefore implies resolving simultaneously for their 
optical properties, in particular the emissivity (k) and its wavelength dependence 
(commonly parametrized as u^). 
Large-scale studies of dust temperature show that the grains in starless cores 
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are colder than in the surrounding lower-density medium, as expected from the at- 
tenuation of the ISRF deep inside the cloud. From ISO FIR (170 and 200 //m) ob- 

servat ions toward the vicinity of a number of dense cores. IWard-Thompson. Andre &i Kirk 

found a prevalence of flat or decreasing temperature gradients with cloud 
temperatures of 15-20 K and core values of 10-12 K. Because of the 1 arcminute 
resolution of these observations (0.04 pc at the distance of Taurus), the core value 
represents only an average temperature over the de i isest gas. Also using low- 
resolu tion ISO data and a similar analysis, Toth et al. ( 20041 ) and ldel Bur go Sz Laureijd 
( 2OO5I ) have derived inward temperature gradients that reach central values of 
~ 12 K for a number of Tauru s cores, while for the denser L183 core a value 
close to 8 K has been measured ( Lehtinen et al. 20031 : Pagani et al. 20031 ). These 
temperature estimates, however, ignore the effects of grain evolution and assume 
a constant dust emissivity, while observations are beginning to show that the 
emissivity does in fact vary. FIR observations of a filament in Taurus by the 
PRONAOS balloon-borne experiment suggest that at the highest densities, the 
population of smallest grains disappears and the submillimeter emissivity in- 
creases by a factor of abou t 3, a s expected by models of dust coagulation into 
fluffy grains (ISteonik et al.l l2003l 'l. These observations also suggest inner dust 



temperatures around 12 K, although the 3-arcminute resolution sets an impor- 
tant limitation. 

The trend of decreasing dust temperature with depth observed in the large 
scale maps is expected to continue inside the dense cores, but the exact shape 
of the resulting gradi ent is still a matter of some unc e rtaint y. Initial models 
by lEvans et al.l (12001^ and IZucconi. Walmslev fc Galiil (|200lh suggested sharp 
core edge-to-center drops in temperature of about a factor of 2 (an example 
is provided in Fig. Uj). However, these models either assumed that cores are 
directly exposed to the ISRF (correct for Bok globules on ly), or use a simple 



attenuation law for cores that are embedded in clouds. lAndre et al. (I2OO3I ) 
and IStamatellos k Whitwor"thl (|2003l ) have argued that the effect of the sur- 
rounding cloud is to flatten the temperature gradient, so that at the center of 
an embedded core the gradie nt is less pronou^lced fsee Istamatellos et a,l. l B: 
Goncalves. Galli &: Walmsleyll2004l . for 2D and 3D radiative transfer models). Ir- 



respective of the models, however, central dust temperatures close to 7-8 K are 
predicted for typical dense cores, with the exact value depending on the dust 
central density and the presence or not of a surrounding cloud. 

From an observational point of view, possible gradi ents in the du s t emi ssivity 
still limit our estimate of the temperature profiles. iKramer et al. ( 20031 ) have 
combined submm emission and NIR extinction data to study several starless 
cores in IC 5146. With an angular resolution of 30 arcseconds, these authors 
find that the dust temperature is typically 12 K and close to constant inside 
the cores, while the surrounding, lower extinction material has a higher value of 
about 20 K. They also find that the dust emissivity increases systematically by a 
factor of 4 between 20 and 12 K, consistent with models of dust evolution (e.g., 
Qssenkopf &: Henningj IT994I ) . Using similar technique and angular resolution. 



Bianchi et al.l (|2003l ) have also found an emissivity gradient in the B68 globule, 



together with an outer d ust temperature of 14 ifc 2 K and an inner temperature 
drop of a factor of 1.5. Schnee Goodman ( 20051 ). on the other hand, have 
studied the submm emission from TMC-IC with 14" resolution and found a 
central gradient in the submm color. As illustrated in Fig. \5\ this color gradient 
can result from either a central temperature drop to about 7 K or a systematic 
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change in the dust emissivity index (or both). Clearly more data and modeling 
are needed to characterize the existing gradients in dust conditions inside dense 
cores. Future multiwavelength mapping in the FIR with the Herschel satellite 
offers the best hope to finally disentangle the changes in temperature and optical 
properties. 

3.2.2 Gas Temperature Like the dust temperature, the temperature of 
the gas is determined by the balance between heating and cooling. For the 
lower density cloud gas, heating mostly occurs through ionization by cosmic rays, 
while the cooling is rnainly due to line radiation from molecules, especially CO 
( Goldsmith &: Lange il ll978l l. At the high densities of cores (more than a few 10 
cm~^), gas-dust coupling via collisions becomes important, and this process will 
additionally heat or cool the gas depending on the difference between the gas and 
dust temperatures ( Burke &: Hollenbach 19831 ). Also at these densities, and at 
typical core temperatures, CO and other molecules start to disappear from the gas 
phase due to their freeze out onto the dust grains (^ 14. 2p . diminishing the ability 
of the gas to cool down via line radi ation. The detailed study of thermal balance 
in dense cores bv iGoldsmithI ^200l\ ) shows that the gas-dust coupling, together 
with the high efficiency of the dust radiative cooling, easily compensates the 
disappearance of molecular coolants, and keeps the temperature of the gas deep 
inside the core at rather low values, comparable to those of the dust component. 
At the lower-density core edges, the gas and dust are thermally decoupled and the 
gas becomes warmer than the dust due to photoelectric heating. The temperature 
difference between the two components depends on the strength of the external 
UV field (see Fig. ED. 

Keeping the gas cold and at a close-to-uniform temperature has important 
implications for the mechanical equilibrium of cores, as the gas componeii t is the 
main contributor to the thermal pressure. iGalli. Walmslev fc GoncalvesI (120021 ) 
have studied the equilibrium structure of a pressure-supported dense core with 
a realistic treatment of the gas and dust cooling. Although the resulting gas 
temperature distribution is not exactly constant, its gradient is shallow enough 
to require for stability a density profile which is very close to the Bonnor-Ebert 
solution (see Fig. |3|). Thus, if cores were supported by thermal pressure only, the 
BE model would be consistent with the internal thermodynamics. 

To measure the temperature of the gas component in clouds and c ores, the leve l 
excitation of simple molecules like CO and NH3 is commonly used (|Evanslll999l ). 
CO and its isotopologues have small dipole moments, so their lowe r energy level s 
thermalize at relatively low densities of the order of few 10^ cm~^ ( Spitzer 19781 ). 
typical of the extended gas in dark clouds. Temperatures derived using CO in re- 
gions with no star formation indicate gas temperatures of 10-1 5 K, with a possible 
increase towards the lower density gas near the cloud edges (ISnel]|ll98ll ). Tem- 
peratures for the denser gas in cores are better determined using the metastable 
inversion levels of ammonia ( NH3), which can be observed simultaneously at a 
wavelength of about 1.3 cm (jWalmslev fc UngerechtsI Il983l ). NH3, in addition, 
appears to rise in abundance in the densest (nHa > 10"^ cm~^) regions (§112]), so 
it tends to selectively trace dense core gas. Low angular resolution NH3 obser- 
vations of dense cores typically indicate gas temperatures of about 10 K with a 
rather narrow range of variation for clouds like Taurus, Perseus, an d Ophiuchus 
(jMyersI Il983l : iBenson Myerd Il989l : Ijijina. Myers fc Adariisl Il999l ) . These low- 
resolution measurements represent core-wide averages of the temperature, and 
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are in reasonable agreement with the theoretical expectations mentioned before. 
Due to the relatively low angular resolution of the NH3 observations (caused by 
the low frequency of the transitions) , few studies of th e radial dependence of the 
gas temperature inside dense cores have been made. iTafall a et al. have 
studied the temperature profiles of two Taurus- Auriga starless cores, L1498 and 
L1517B, using NH3 observations with a resolution of 40" (0.027 pc). In both 
cores, the temperature remains constant over the central 0.1 pc with a value of 
about 10 K. Higher resolution NII3 observations of the denser starless core L1544 
seem to indicate a decrease of the gas temperature toward the innermost parts of 
the core (Cr apsi et al. in prepar ation). For the starless Bok globule B68, on the 
other hand, iBergin et al.l (|2006l l have found an opposite effect: while the inner 
core traced by NH3 is at about 10 K, the surrounding layers, traced with CO 
isotopologues, are slightly cooler (7-8 K). This central warming may result from 
a weakening of the gas-dust coupling at the highest densities due to grain coag- 
ulation. Given this diversity of results, further observations of the temperature 
gradients inside dense cores are needed to clarify the origin of these significant, 
but relatively small (20%), temperature differences. 



3.3 Velocity Structure 

In contrast with the supersonic velocity fields characteristic of the clouds as a 
whole, dense cores have low velocity internal motions. Starless cores in clouds like 
Taurus, Perseus, and Ophiuchus systematically present spectra of core-tracing 
species (e.g., NH3) tha t have close-to-thermal linewidths, even w hen observed at 
low angular resolution ( MversI 19831 : Jiiina. Mvers &: Adamsiri999l ). This indicates 
that the gas motions inside the cores are subsonic, either in the form of random 
turbulence or as part of more systematic global patterns (Fig. [6|) . 

Despite the slow speeds involved, the internal motions of cores are in general 
complex, with no single element identified as dominant. The excess in linewidth 
over the thermal value has commonly been interpreted as resulting from turbulent 
motions, although as we will see below, other velocity patterns like infall could 
contribute to the line broadening. Even in this interpretation, turbulent motions, 
by being subsonic, contribute less to the gas pressure than the the rmal compo- 
nent, so they represe r it a minor contribu t ion to the core support ( MversI 1983 ; 
Goodman etal]ll998l : ICaselh et ai]l2002al : iLada et al.ll2003l : iTafalla et al.l 120041') 



Studying the spatial distribution of turbulent motions in cores, iGoodman et al. 



fl998) have fo und that they deviate from the linewidth-size relation seen at large 
scales (§ ESI lLarsonllT98lh . Instead, the non-thermal component of the NH3 
linewidth in a c ore is approximately constant inside a " coherence" radius of about 
0.1 pc (see also Caselli et al.ll2002a : Tafalla et al.l 2004 ). This coherence radius is 
close to the 0.04 pc radius found to separate the binary and clustering regimes in 
Taurus and interprete d as the lengt h-scale at which the cloud self-similar behav- 
ior breaks down (^ 12. 5^ Larson 19951 ). It is therefore tempting to suggest that the 
coherence radius is the scale at which turbulence loses its grip on cloud motions 
and allows the formation of thermally dominated cores. 

A rotation components is expected in any objec t that forms by the gravita- 
tional contraction of galactic material ( Spitzer 19781 ). From the study of the spa- 
tial changes in the velocity centroid of lines, however, it has been known for some 
time that dark clouds as a whole are not rotating rapidly ( Heiles &: Katj 1976 : 
Arauilla fc Goldsmithlll986l l . Observations of dense cores and globules reveal a 
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similar lack of fast rotation in the dense gas. iGoodman et al. (|l993l ) analyzed the 
rotation patterns of 43 dense cores as derived from observations of the NH3(1,1) 
line, and estimated that the typical ratio between rotational and gravitational 
energies is only 0.02. These authors also found no correlation between the ro- 
tational axis and the orientation of core elongation, which is again indicative of 
rotation being only a minor element in core suppor t. From higher reso lution 
N2H"'' observations of an equivalent sample of cores, ICaselli et al. have 
derived similarly low values of the rotational energy, and found that the point-to- 
point velocity gradients inside the cores are much more complex than expected 
for simple solid-body rotation. Rotation is therefore only one of the ingredients 
of core kinematics. 



The slow rotation rate of dense c ores is in fact remarkable. lArquilla Sz Goldsmith 
(|l986l l and lGoodman et al.l (|l993l ) have shown that the specific angular momen- 
tum j in a core or a cloud scales with its radius R approximately as R^'^. Clouds 
of 1 pc radius have typical j values around 10^^ cm^ s~^, while typical cores 
(0.05 pc radius) have j values of about 10^^ cm^ s~^. If this trend of decreasing 
j reflects the evolution of specific angular momentum as the cloud gas undergoes 
contraction, an efficient mechanism of angular momentum removal is required 
to reproduce the observations, and in particular, to explain the low j values of 
cores. Magnetic coupling between the contracting gas and the lower density en- 
vironment has been classically favored as the mechanism for transferring angular 
momentum from a contracting c ore to its surroundings ( Mestel &: Spitzeil 1956 : 
Mouschovias &: Paleoloeo ^ Il979l ^. although an alternative int erpretation in terms 



of gr avoturbulent fragmentation has recently been proposed (jJappsen &: Klessen 
2OO4I ) . The low j value of cores is in fact comparable to or slightly larger than 
the value inferred for protostellar envelopes, T Tauri disks, and wide binaries, 
suggesting that the later phase of gravit ational contraction of the core occurs 
with almost constant angular momentum ( Bodenheimer 19951 ). Of course, a final 
phase of angular momentum redistribution via disk transport or tidal interactions 
between stars is still needed to e xplain the low specific angul a ,r momentum found 
in individual T Tauri stars (e.g., Shu. Adams &: Lizanol 1987 : Bodenheimer 1995 : 
Larsonl[200i ). 

Inward motions can also contribute to the dense core kinematics. They are 
inferred from the observation of optically thick, self-absorbed lines of species like 
CS, II2CO, or HCO'^, in which low-excitation foreground gas absorbs part of the 
background emission. In a core with inward motions, the foreground gas moves 
toward the background, so the absorbing material is red-shifted with respect 
to the emitting layers. The resulting spectrum is a self-absorbed profile with 
a red-shifted dip and a relatively brighter b lue peak (see CS J=3-2 in Fig. [6l) 



a pattern often called an infall asymmetry (jLeung Sz Brownl 119771 : 1 M vers et al. 

In a starless core, this spectral feature has been first identified towards 
L1544, for which a.n average inward velocity of about 0.1 km s~^ was inferre d 
(|Mvers et al.lll996l : iTafalla et all 19981 : IWilhams et al.lll999l:[Caseni et al.ll2002bl') 
A foll ow up single-pointing survey of 220 starless cores by iLee. Myers &: Tafalla 
(|l999l l has found a statistically significant over-abundance of inward motions 
with typical velocities of 0.05-0.1 km s~^, indicating that subsonic contraction is 
a common feature of core kinematics (also lGregersen fc Evand[2OO0l ). The spatial 
distribution of these motions is typically large, 0.05-0.15 pc, and comparable t o 
the observed size of the cores (|Tafalla et al.l Il998l : iLee. Mvers fc Tafallal l200lh . 
Such a large large extent could result either from the gravitational collapse of a 
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core having a centrally condensed ini tial state ( Mverd 20051) or be a ssociated with 
the process of core condensation (see iMvers. Evans &: Ohashill2000l . for a review) . 

Additional velocity patterns ha ve been proposed t o explain the observations of 
molecular lines in different cores. iLada et "ah have found an asymmetric 

velocity field with both inward and outward motions in the B68 globule, and they 
have interpreted it as resulting from small amplitude pulsations of the core outer 
layers. Such a pattern seems in fact co nsistent with the cl ose to Bonno r- Eber t 
density structure found in this object ( Keto et"al] 20061 ). Tafalla et al. ( 20041 ) 
have proposed that residual internal motions in L1498 and L1517B may originate 
from asymmetric contraction, as the motions are correlated with asymmetries 
in the pattern of depletion of species like CO and CS. Finally, by comparing 
the velocity centroids of NgH" *" and CO isotopologues in a sample of 42 cores, 
Walsh. Mvers &: BurtonI ( 20041 ) have concluded that the relative velocity between 
the dense cores and the ambient medium is low (< 0.1 km s~^). This implies that 
newly born stars are not expected to move appreciably from their natal cores. 



3.4 Magnetic Field 

Estimates of the magnetic field in starless cores are mainly obtained from Zeeman 
splitting measurements and from polarization observations of the mm/submm 
dust continuum. Zeeman measurements have the disadvantage of a low angular 
resolution and the use of tracers sensitive to molecular depletion, so they are prone 
to miss the densest part of the cores. On the other hand, they are unique because 
they provide a direct estimate of the strength of the (line of sight component of 
the) magnetic field, -B/ n,;. From OH Zeemari measurements of the L1544 core with 
the Arecibo telescope, Crutcher fc Trolandl (|2000; ) have estimated a line of sight 
strength of about 11 //G averaged over a beam of 3 arcmin (about 0.1 pc), in 
good agreement with the prediction from a previous model of ambipolar diffusion 
(ICiolek fc B"^l200nl ). These observations, however, are more sensitive to the 
large-scale cloud gas than to the core material because of the the low dipole 
moment and the chemical distribution of the OH molecule in clouds (§ 14. 2p . 
Zeeman measurements using tracers of dense gas have mostly produ ced negative 
results, like Bios < WO uG in L1498 (CCS data, iLevin et allboOlh and Bin. < 
15uG in TMC-1 (C4H data. lTurner fc Heilegll200fil ;i. IShinnaga. Tsuboi fc Kasuga 
(I1999I ). on the other hand, have reported an estimate of 160ib42;uG in the young 
starless core L1521E from CCS observations. 

Submm polarization measurements provide a direct estimate of the orientation 
of the (plane of the sky component of the) magnetic field. They can also be used to 
estimate indirectly the strength of the field by cor nparing the spatial dispersion o f 
its orientation and the turbulence level of the gas ( Chandrasekhar Sz Fermi 19531 ). 
As mentioned in § 12.71 these measurements often lose sensitivity towards the 
brighte st, highest density regions, where the polarizat i on fraction drops with in- 
tensi ty (|Matthews fc Wilsonll200ol : iHenning et al.ll200ll : IWolf. Launhardt &: Hennine 
20031 ). which again makes the estimates more sensitive to the field in the outer 
parts of the cores. Because of the combination of weak submm emission of the 
starless cores and a low polarization fraction (~ 10%), these observations are 
challenging and have only be en carried out in a redu c ed nu mber of objects. For 
cores L1544, L183, and L43, IWard-Thompson et al.l (|2000l ) and ICrutcher et al. 



( 2OO4I ) have found polarization vectors that are fairly uniform but not aligned 
with the minor axis of the cores by about 30°, as it would have been expected 



Cold Dark Clouds 



25 



by ambipolar diffusion models (but see Basu 2000l for an alternative view if the 
cores are triaxial). From the dispersion of the polarization vectors and the turbu- 
lence level of the gas, these authors estimate plane-of-the-sky intensities between 
80 and 160 /iG. Strictly speaking, these field intensities make the cores super- 
critical (i.e., with a magnetic field too weak to balance gravity). However, the 
inclusion of a geometrical correction factor moves the cores to the critical or 
slight ly subcritical regime , suggesting that the field may be dynamically impor- 
tant (ICrutcher et al.l 120041 ). Lower field values of 10 and 30 fiG have been esti- 
mated towards the cores L1498 and L1517B, respectively, indicating that these 
cores may be supercritical by factors of ab out 2 even after geometrical correction 
( Kirk. Ward-Thompson Crutched 120061 ). 

As with the large scale magnetic field measurements (§ 12. 7p . the observations 
of cores provide an ambiguous answer to the question of whether magnetic fields 
dominate or not the equilibrium of dense cores. Clearly more observations of both 
the Zeeman effect and the dust polarization are needed to clarify the role of mag- 
netic fields in the process of core and star formation. The close to critical values 
observed so far, however, suggest that unless observations have been systemati- 
cally biased in one or the other direction, cores appear to be close to the point 
of criticality to within a factor of about 2. Explaining this apparent coincidence 
should therefore be a main requirement to any theory of core formation. 



4 CHEMISTRY OF PRE-STELLAR MOLECULAR CORES 



4.1 Background 

The presence of gas phase chemical inhomogeneities within molecular clouds 
has been known for some time. In particular, spatial differences in emission 
morphologies between carbon chains (e.g. HC3N, C2S) and nitro gen- hydrogen 
mole c ules (NH3, N^H+) on scales of 0.05-0.2 pc were described (jSuzuki et al. 



I992I : IPratap et al1ll997l l. Surveys of numerous cores in Taurus revealed larger 
emission sizes and vel ocity linewidths for CS when compared to NH3, despite 
similar line centroids ( Zhou et al. 19891 ). These differences have generally been 
attributed to the slower timescales to activate the nitrogen chemistry as opposed 
to the carbon chemistry, or perhaps due to sor ne dynamical cycling of mate- 
rial that can enha nce the carbon chemistry (see Ivan Dishoeck and Blakd Il998l : 
Langer et al.ll2000l . for a discussion of general chemistry). 

In this regard two low-mass cores, TMC-1 and L134N, were isolated and stud- 
ied as template objects because of their rich chemis try and due to the fact that 



they are unassociated w ith star formation activity (jSwadd Il989l : iHirahara et al. 



I992I : IPratap et al.lll997l ). Over one thousand papers have been written for these 
two objects alone. In general chemical s tudies comp u ted a verage line of sight 
abundances using a technique outlined in llrvine et al.l (|l987l V H2 does not emit 
appreciably at the cold temperatures associated with these objects (T ~ 10 K). 
Therefore optically thin isotopologues of CO, generally C^^O, were used as sur- 
rogates to trace the H2 column density with previously calibrate d abundances 
(e.g. x{C^^O) ~ 1.7 X 10"'^ = nfC^^0\/nfH9l: lFrerking et al.lll982l '). These stud- 
ies strongly aided the interpretation that ch emical reacti ons between ions and 
molecules dominate the gas-phase chemistry ( Herbst 20051 ). 

In parallel to gas-phase chemical studies, the fundamental vibrat ional modes of 
solid state molecules such as H2O and CO were detected in the ISM ( Gillett &: Forrest 
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19731 : iLacv et al.lll984l ). Many studies of molecular ices frozen on the surfaces of 
dust grains have been conducted towards embedded sour ces which provide a 
strong infrared background to study foreground ices (e.g. Ivan Dishoeckl[2003 ). 
To provide important insight on ices in cold gas, observations of bright field stars 
located behind molecular cloud material (with extinctions < 20 mag) are used a s 
candles that probe material remote from embedded sources ( Whittet et al. 19881 ). 
Icy mantles coating grains in ambient g dominated by H2O ice with sub- 

stantial amo unts of CO and C O? at a level of ~25% of the water ice abundance 
in each case ( Gibb et 311120041 ). 



There exists a threshold extinction below which the ice features are not seen, 
implying that grains are not mantled at low extinction. This threshold varies 
for each ice component and is different between various clouds. In Taurus, the 
water ice threshold is Ay(H20) = 3.2"* it 0.1"*; a compar able threshold is ob 



served for CO9. ice but for CO A^^(C O) = 6.8™ ± 1.6"* (jWhittet et al.l I2OOI 



Bergin et al. 20051 : Whittet et al.|[2007l ). The threshold value presumably is re 



lated to the formation mechanisms for ices and in Taurus the relatively low water 
ice threshold may imply that the water ice forms at a low density were atomic 
H is available. These observations further demonstra te that ices constitute a 
substantial fraction of the available carbon and oxygen (jvan Dishoeck and Blakd 
[l998). Thus CO (and its isotopologues) is not a reliable tracer of the H2 column 
when using previously calibrated abundances. In addition, it is clear that grains 
can act as a catalyst for chemistry on the surface, producing both simple and 
complex molecules as part of an interchan ge between the gas and solid phases 
(iTielens fc Hagenlll982l : iHerbst et alll200,5l ). 



4.2 Observational Chemistry 

4.2.1 Gas-Phase Freeze-Out As discussed in §3 a sample of low-mass 
starless cores are now known with well described physical structure as a function 
of position: Njf2(r), T^Hiir), Tdustir), Tgas{r). Sample profiles are given in Fig.H 
Foremost among these are the two new template objects: L1544 and B68. There 
are some caveats to the determination of each of these parameters, but the overall 
structure is well characterized and has led to significant advances in our under- 
standing of the chemistry. In general there are two factors that influence these 
gains. (1) Using the dust mm/submm continuum emission as a surrogate for the 
II2 distribution allows for an examination of molecular abunda nce variations di- 
rectly relative to II2. This is a fact noted by lMezger et al. I (I1992I ) in the high-mass 
NGC2024 star forming region. However, the presence of embedded sources can 
substantially inf luence the dust emissioii morp hology, potentially leading to differ- 
ent conclusions ( Chandler Sz Carlstromlll996l ). While starless cores do have small 



thermal gradients (§3.2) this particular complication is minimized. (2) Molecular 
emission depends on the II2 density, gas temperature, and on the abundance. 
Using sophisticated one-dimen sional radiat i on transfer codes (Monte -Carlo or 
Accelerated Lambda Iteration; lBernesl[l979l : iRvbicki k Hummed [l99ll ) and the 
density and temperature structure as inputs, the abundance profile along the line 
of sight can be constrained for the first time. 

Initial studies found the first direct link between chemistry in the gas and 
interactions with grain surfaces by simply comparing the estir nated CO and CCS 
column densities to the H2 column estimated from the dus t ( Kuiper et al. 19961 : 
Willacv et al]ll998l : iKramer et allll999l : ICaselh et aDll999l ). This is interpreted 
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as evidence from gas-phase data that the CO molecules are freezing onto grain 
surfaces. A comparativ e study of CO suggest freeze-out dominates when densities 



exceed ~ 3 x 10^ cm ^( Bacmann et al. 20021 ) 



Spatial mapping of ices within individual cores supports the interpretation of 
the gas-phase observations as the abundance of CO ice is found to significantly 
increase when densities exceed ~ 10 cm" (|Pontoppidanll200fil ;i. Near the core 
center the abundance is close to the typical gaseous CO abundance of 10~^, 
indicating that the majority of CO is frozen out. Furthermore, the abundance of 
water ice is typically 0.5 — 0.9 x 10~^ except at the highest densities where the 
abundance increases ( Pontoppidan et al. 20051 ). The molecular ices are therefore 
substantial reservoirs of the available oxygen. 

4.2.2 Selective Freeze-Out Molecular surveys reveal further differences, 
Fig. [7] shows a sample of this work with a comparison of the 850^m dust contin- 
uum distribution in the Barnard 68 starless core with molecular emission traced 
by C^'^O and N2H"'". In this core the N2H"^ emission more closely follows that of 
the dust than does C^^O, which appea rs as a "semi-ring-l i ke" structure a round 
the dust and N2H"'" emission maximum ( Bergin et al. 2002 ). Tafalla et al. ( 20021 ) 
surveyed 5 cores in Taurus (including LI 544) illustrating general characteristics 
of core chemistry: carbon-bearing species, represented by CO and CS deplete 
from the gas while nitro gen-hydrogen beari ng molecules, N2H"'" and NH3 trace 



the core center (see also lCaselh et al.l Il999l ) . From the abundance profile it is 



estimated that the CO and CS abundance traces a large dynamical range with 
declines of at least 1-2 orders of magnitude in the core centers relative to the 
lower density core edge, while the abundances of nitrogen molecules either stay 
constant or decay more slowly. The interpretation of "selective" freeze-out, where 
molecules exhibit different behavior in terms the response of the chemistry to in- 
teractions with grain surfaces, naturally explains some long-standing issues such 
as the emission di fferences seen between carbon an d nitrogen molecules described 
earlier (compare ISwadel 19891 : Pagani et al.ll2OO.5l ). and the larger core sizes and 
velocity dispersio i is seen for CO an d CS when compared to earlier NH3 data 
(jZhou et al.lll989l : iTafalla et all I2OO2I ) . For star formation this isolates the nitro- 
gen hydrides as the key probes of dense gas. 

4.2.3 Carbon, Oxygen, AND Nitrogen Cores with well defined physical 
properties have clarified additional outstanding chemical issues regarding carbon, 
oxygen, and nitrogen chemistry. It is clear that CO represents the major reser 



voir o f carbon in the gas and, along with CO2, in the grain mantle. iMaret et al. 



(|2006l l demonstrated that N2 is not the dominant carrier of nitrogen in molecu- 
lar gas. They suggest that N I as a main reservoir of atomic nitrogen (see also 
van Dishoeck and Blakd 19981 ): however, at present, it is unclear what is the pre- 



dominant form of gas-phase nitrogen. Observations by the Submillimeter Wave 
Astronomy Satellite (SWAS) and Odin also demonstrated that the abundance of 
ortho-water vapor and n iolecular oxygen i s well below theoretica l expectations 
throughout the Galaxy (ISnell et al.1 l20qol: iHialmarson et aD I2OO5I) and in star- 
less cores in particular (jBergin fc Snel]ll200l l~ Due to atmospheric absorption, 
para-water vapor has yet to be observed and the limits on the water abundance 
disagree with pure gas-phase chemical theory provided that the ortho/para ratio 
is not below the equilibrium value at 10 K (0-H2O/P-H2O ~ 0.3). The low water 
vapor abundance is interpreted as the result of water ice formation on grain sur- 
faces. At grain temperatures below ~ 110 K (§4.3.2) water will remain frozen on 
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grains, which indirectly lowers the gas phase atomic oxygen abundance thereby 
hindering formation o f water and molecular oxygen in the gas (jBergin et al.ll200d : 



Charnlev et al. 200ll ). Further analysis suggests that the water vapor emission 



likely arises from the surface photodissociation region where a balance between 
photodesorption of water ice and p hotodissociation of wat er vapor exists (OH 
should behave in a similar fashion) ( Melnick &: Berg"inl 20051 ). The launch of the 
Herschel Space Observatory, which will be capable of observing both ortho- and 
para-water will provide definitive information on the water abundance and its 
potentially defining role in gas-phase freeze-out. 



4.3 Chemical Theory 

4.3.1 Basic Processes In the literature the loss of gaseous molecules to 
the solid phase is often labeled as depletion, which implies a lowering of the gas- 
phase abundance, or as freeze-out, which is more direct statement of the process. 
In practice both terms can be used, but there are some nuances. While neutral 
species freeze onto grain surfaces it is not thought that molecular ions behave in 
a similar fashion. Grains in dense clouds are believed to carry negative charge 
( Weingartner &: Draind boOll ). The ion- negatively charged grain recombination 
releases a few eV which, if carried by the prod ucts, significantly exceeds th e 
molecule-grain binding energies of ~ 0.1 - 0.5 eV (IColhngs k McCoustral[2005l ). 
Thus ions can deplete, via, for example, gas phase processes and by pre-cursor 
molecules (e.g. CO for HCO+) adsorbing onto grains, but do not freeze-out. 

Atoms and molecules freeze onto grain surfaces or adsorb with a rate oi kfo = 
avUgS s~^. Here a is the grain cross-section (generally assuming to be an av- 
erage grain with a 0.1 /um radius), v is the mean velocity of the Maxwellian 
distribution of gaseous particles, Ug the space density of grains, and S is the 
sticking coefficient (i.e. how often a species will remain on the grain upon im- 
pact). Using relevant parameters for CO at 10 K (assu ming 5 = 1 as suggested 



by laboratory experiments; e.g. Bisschop et al. 20061 ) the freeze-out timescale 



Tfo = kj^ ~ 5 X 10 /nH2{cm~ ) yrs. Once on the grain, molecules can react with 
any species that is mobile on the grain surface at 10 K, such as atomic hydrogen, 
or if a species is saturated (e.g. H2O, NH3, CO2) it will remain inert on the 
surface. For cold dense gas that is shielded from external radiation the physical 
processes that can remove or desorb molecules fro m grain surface are either inac- 
tive or can only desorb the most volatile ices (see Ivan Dishoeck and Blakelll998l . 
for a complete summary of desorption processes). In cold cores it is thought 
that energetic cosmic ray impa c ts can desorb the most volatile ices (e.g. CO 
but not H2O; Leger et al. 19851 ; Hasegawa and Herbst 19931 : Bringa &: Johnson 
20041 ^. Thus the dominant effect is the adsorption of atoms and molecules onto 
grain surfaces and any subsequent consequences that result from this interaction. 
4.3.2 Chemical and Dynamical Timescales In Fig. [8]we compare rel- 
evant chemical and dynamical timescales. Representative dynamical timescales 
are provided for free-fall and ambipolar diffusion , the latter having been su, 



gested as a primary mode of core condensation ( Shu. Adams Sz Lizanol 198?^ . 



The free-fall timescale is rjj = [37r/32Gp]'''^ ~ 4 x lO^/n^^ Y^^ the ambipo- 
lar diffusion timescale is discussed in § I4.4.2[ The gas-phase chemical timescale 
is set by cosmic rays and has a constant dependence on the density, while freeze- 
out has linear dependenc e on the dens ity. The evaporation rate is given by the 
Polanyi-Wigner relation (lTielensll2005l ^: 
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T, 



evap 



(1) 



where is the vibrational frequency of the molecule in the potential well (typ- 
ically ^ 10^^ s^^) and Ef, is the binding energy. The evaporation timescale has 
an exponential dependence on the dust temperature, but critically depends on 
the binding strength to the grain surface or man tle. Measured binding ener- 
gies for CO-CO ar e ~ 820 K lCoUings et al.ll2003l l and ~ 5800 K for H2O-H2O 
(jFraser et al.l 120011 ). implying evaporation temperatures of ~ 15 K for CO and 
~ 110 K for H2O at 11^2 = 10^ cm~'^ and assuming a typical 0.1 fim grain abun- 
dance of Xgr = 10"^^. Fig. [8] demonstrates a key aspect of starless cores - as cold 
objects, T^ust ~ Tgas ~ 10 K, evaporation is unimportant, with timescales > 10^ 
yrs. Furthermore the freeze-out timescale becomes shorter than the free-fall and 
ambipolar diffusion timescales as the density increases. Thus as the core physi- 
cally evolves, the effects of freeze-out on gaseous species are magnified, and can 
be used to constrain the mechanism of core formation. 

4.3.3 Chemical/Dynamical Models Combined chemical and dynamical 
models follow the chemical composition as the core contracts. Such models gen- 
erally adopt parameterized fits or a more direct coupling to simple dynamica l 
models (e.g. iLarsonl [l969l : IPenstonI [l 969|: IShulll977l: iBasu MouschoviasHigoJ) 



Notab le chemical /dynamica l mod els are [Rawlings et al.l (|l992l ). iBergin &: Langed 



(119971 ). lAikawa et al. 



^ (|200lh and iLi et al. (2002). With bas ic gas- gram micro- 
physics incorporated into time-dependent models fsee lTielenslboO.^I ;) the chemical 
structure discussed in § l4.2.2l can be explained and, in some cases, was predicted. 

Basic elements of the chemistry are illustrated in Fig. [9l During collapse or 
condensation, the central regions of the core increase in density to levels where 
neutrals rapidly freeze onto grains (n^/j > 10^ cm~^). This contrasts with the 
lower density less evolved core edges with undepleted abundances. As collapse 
proceeds t hese effects becoine amp lified (compare the two times provided for CO 
in Fig. [9|). Bergin &: Langei] (1997) first showed that carbon- and sulphur-bearing 
species (e.g. CO, CS, CCS) suffer the greatest effect fro m depletion. In con trast 
the nitrogen hydride pool is the least affected (see also Aikawa et al.ll200"ll ). At 
first this was assumed to be due to a lower binding strength of N2 (the precursor 
molecule for both N2H+ and NH3) to the g rain surface than for CO. This has 
now been measured to be equivalent to CO (j Oberg et al. 20051 ). 



It is now recognized that the observational appearance of selective freeze-out 
can be separated into first-order chemical effects, the direct freeze-out of neutrals 
(e.g. CO, CCS, CS, N2), and second-order effects, changes created by the im- 
balance in the chemistry left primarily by the loss of gaseous carbon monoxide. 
CO is a major destroyer of molecular ions and its removal fror n the gas leads 



to a change in the relative abundance of major charge carriers (| Rawlings et al. 
1993). Particularly affected is the Kg ion, the precursor to molecular ions such 



as N2II"'" and HCO"*" that are the respective daughter products of N2 and CO. 
These key observable ions (N2H"'" and HCO"*") will also decrease in abundance as 
the parent molecules freeze-out; however, CO is a major destroyer of N2II"'" and 
as CO disap pears from th e gas there is a subsequent increase in a. bundance of 



N2H+ (, Aikawa et al.ll200ll: iBerdn et al.l[200i: I.Torgensen et al.l l2004l'l. NH3 then 



forms from NoH+ (jCeppert et alJl2004l : lAikawa et al.ll200,5l ). Thus the abundance 
of the nitrogen hydrides is strongly enhanced and these species are probing the 
gas where CO (and other carbon-bearing species) are depleting. An important 
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point here is that the abundance of both NH3 and N2H+ are several orders of 
magnitude below the abundance of molecular nitrogen. N2 is freezing onto grains 
at the same level as CO. Provided the N2 abundance does not deplete to levels 
below N2H"'" and ammonia, then these tracers can still form and emit from lay- 
ers where commonly observed tracers of dense regions have significantly reduced 
abundances. 

Combined models reproduce not only the pattern of depletions, bu t also the 
column densities observed in L1544 (lAikawa et al.1 120011 : iLi et al.ll2n02l ;i. In gen- 
eral models obtain best agreement with collapse timescales much shorter than 
would pro ceed from ambipolar d iffusion starting from a magnetically subcritical 
state (e.g. IWalmslev erall l2004l l : however there is some disagreement (jLi et al. 



2002! ). A clear point is that evolutionary timescales at densities near ~ 10' 



cm 



cannot be much longer than ~ 0.5 — 1 Myr as si gnificant CO freeze - out would be 



prod uced at densities below where it is observed (jFlower et al.ll2005l : iTafalla et al 




Moreover, collapse cannot be too fast (> 0.2 km s ^ at Tdyn ;^ 5 x 10^ 



o r ther e will be less depletion and emission linewidths would be too broad 
3[l98i). 

Current work has greatly benefited from laboratory work on lo w temperature 

gas phase reactions, grain binding energies, and on surface catalysis ( Florescu-Mitchell Sz Mitchelll 
2006; Collinas & McCoustra 2005: lNagaoka et al.l 2005: Bi sschop et al.ll2006l : lHornekaer et al. 



2OO3I : iKatz et al.lll999l ). Moreover the recent re-calculation and extension of col- 



lisional r ates of astrophysically impo r tant molecules is of great use a, nd deserves 
mention (iDubernet &: Grosjean I2OO2I : lFloweill200il : IPaniel et ai]l2005l ) . However, 
large questions remain regarding the nature of the grain surface catalysis in this 
cold environme nt and the strength of the mechanisrns that could return mol ecules 
to the gas (e.g. iLeger et al.lll985l : iBringa &: ,Tohnsonl[2003 : IShen et al.l[200i ) . It is 
also unclear how complex molecules form in the gas-phase. For example CH3OII 
is detected in cold clouds, but currently has no known gas-phase formation route. 
An examination of the extent of the complex chemistry in cores with well defined 
properties is needed and will likely lead to gain s in our understanding of inter- 
stellar organic chemistry (see Hirota et al.l 2004 ) . 

One clear positive statement is our knowledge of the chemistry and the time 
sensitive nature of the gas-grain interaction is such that we n ow have the meth- 
ods to place relative ages between cores within a given sample ( Hirota et al.|[2002 ; 
Tafalla &: Santiago 20041 ) . This should aid in our analysis of the core/star for- 



mation process. More over, the combi nation of chemical/dynamical models with 
radiative transfer (see lLee et al. 2004, as an example) offers the best opportunity 
to take these models to the next level and truly begin to test star formation 
theory. 



4.4 Deuterium Chemistry and the Ionization Fraction 

4.4.1 Deuterium Fractionation Another second-order effect induced by 
the CO freeze-out is a sharp rise in the level of deuterium fractionation. Enhance- 
ments of deuterium bearing molecules relative t o the hydrogen counterparts has 
been known for some time (iGuelin et al.lll982l ). In cold clouds enrichments of 



^The typical non-thermal linewidth in starless cores (e.g. L1498, L1517B, B68, ...) is 0.2 
km/s. If this represents a purely inward motion, it would correspond to an infall speed of 0.1 
km/s at most. With a typical core radius is 0.05 pc, this implies a typical contraction time of 
0.5 Myr. 
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2-3 orders of magnitude are observed above the atomic hydr ogen value of (D/H ) 



> (2.3 ± 0.2) X 10-5 estimated within 1 kpc of the Sun (iLinskv et al.ll2006l ). 



Because of the lower zero point energy of deuterium bonds compared to bonds 
with hydrogen, ion-molecule reactions in the dense ISM are th ought to be the 
mechanism responsible for these enrichments ( Millar et al. 19891 ). Deuterium en 



richments can also be created via reactions on the surfaces of dust grains (iTielensI 
19831 : iNagaoka et al.l[2005l ;i. but in these cold cores freeze-out is believed to dom- 
inate as opposed to sublimation (e.g. freeze-out dominates the chemistry of the 
volatile CO molecule). 

In the cold dense ISM, deuterium chemistry is driven by the following reaction: 



H+ HD ^ H2D+ + R2 + 230K. 



(2) 



The forward reaction is slightly exothermic favoring the production of II2D+ at 
10 K, enriching the [D] / [H] ratio in the spe cies that lie at the heart of interstellar 
ion-molecule chemistry ( Millar et al. 19891 1^. These enrichments are then passed 
down the reaction chains to species such as DCO"*", DCN, HDCO, and others. 

Pure gas-phase mo dels without freeze -out cannot produce significant quanti- 
ties of doubly (NIID2 ; iRouefF et al. 2000l ) and triply deuterated ammonia (ND3 ; 



Lis et al.l[2003 ) as observed in starless cores. This motivated a re-examination 



of the basic deuterium chemistry. The primary advance in our understanding is 
two-fold: (1) deuteration reactions do not stop with II2D"*", rather they continue 
towards the formation of both D^H^ and Pj", via a simila.r reac tion sequence 



(IPhilhps fc Vaste]ll2003l : Roberts et al.ll2003l : IWalmslev et~aDl2004l ;i: 



H2D+ + HD ^ D2H+ + R2 + 180K, 



(3) 



H2D+ HD ^ Dt -h H2 + 230K. 



(4) 



(2) The freeze-out of CO, which is primary destroyer of both Hg" and 1120"*", 
increases the rate of the gas phase fractionation reactions (see deute rium species 
in Fig. [9J IStark et al.lll999l : lAikawa et aTll200ll : iBacmann etalll2003l ) . This infer- 
ence is strongly supported by the detection of H2D^ and D2H ^ in starless cores 
(jCaselh et alj|2003l : FVastel et al.|[20o3 : Ivan der Tak et al.ll2005l ). as illustrated in 
Fig.ini 

The creation of multiply deuterated ions leads to the production of high deuter- 
ation levels similar to those observed. An additional by-product is an enhanced 
atomic (D/H) ratio which will accrete onto gra i ns and be available for sur- 
face reactions (lTielensl[l983l : ICharnlev et al.lll997l : IStantcheva and HerbstI l2003l : 
Aikawa et al.ll2005l ) creating molecules rich in deuterium to be exposed when ices 
sublimate off warm grains as the result of stellar birth. However, the extent of 
the surface deuteration opposed to the gas phase needs to be further explored. 
It is worth emphasizing that the observed form of H2D+ in interstellar clouds is 
ortho-H2D+ and models show that the H2D"'" ortho/para ratio (o/p) will change 
with density and is a se nsitive function of the molecular hydrogen ortho/para 
ratio (jFlower et al.ll2006l ). If the H2 o/p ratio is not in equilibrium at the low 
Tg < 20 K (wh ere o/p <^ 1) and instead is near 3:1, then deuteriu m fractionation 
will not work ( Pineau des Forets et al.lll991 : Gerlich et al.l 2002 ). Observations 



■^Reaction [2] has been measured in the lab at lo w temperatures fin ding that there is an 
additional dependence on the ortho/para ratio of H2 (|Gerlich et al.ll2002l ). 
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from the Spitzer Space Telescope suggest t hat the H2 o/p ratio is in equihbrium 
in cold pre-shock gas (jNeufeld et alj l2006l ). a fact consistent with the observed 
fractionation. Herschel and SOFIA will be needed to constrain the H2 o/p ratio 
via observations of ortho- and para-H2D"'" in starless cores. 
4.4.2 Ionization Fraction For subcritical starless cores (where the mag- 
netic field supports the cloud against gravity, 12. 7p the evolutionary timescale 
is set by ambipolar diffusion and t he ion f ractio n is a critical factor in deter- 
mining the magnetic flux leakage. iMcKee (|l989l l discussed the various charge 
carriers and recombination rates deriving a basic expression that has wide use: 

Xp = 1.3 X 10~^?T.(H2)~'^'^. The effects of freeze-out can alter both the pre-factor 

J 11 — — ' ' — — ■ ' 

and the exponent of the dens ity dependenc e (iBergin &: Langerl l997: C aselli et al. 

2002d : IWalmslev et al.ll2004l ). but see also ICiolek Mouschovias, (JmS)- 



Past studies of the ion fraction in dense cores focused on the DCO^/HCO"^ 
ratio, that was predicted to be a sensitive function of the electron abundance. 
Furthermore, it was thought that HCO^ was a m ajor charge carrier throughout 
the core (e.g. ICaselU et al]ll998l : IWilhams etaP [19981 . These studies found line 
of sight average ion fractions of Xe ~ 10~^ (relative to H2), implying weak cou- 
pling to the magnetic field. One variable assumption in models is the presence 
or absence of metal ions (e.g. Fe"*", Mg"*", ...), which have slower recombination 
timescales than molecular ions. If present, metals could dominate the electron 
fraction over molecular ions ultimately produced by cosmic ray impacts. Cur- 
rently detailed models i mply no significant contribution from heavy metal ions 
to the electron fraction flCaselh et al.ll2002cl : iMaret et al.ll2006l ). 

Due to the freeze-out of neutrals, HCO"'" does not trace the core center and 
major charge carriers in the core center are likely to be a combination of Hg", 
1120"*", D2H'*', and Dg". An example of the "typical" ionization structure is given 
in Fig. [TTJ When illuminated by the standard interstellar radiation field or below, 
the ionization is high at core edges {Ay < 1-2 mag) due to CO photodissociation 
and ionization of carbon. Deeper in the cloud the regime of cosmic ray ionization 
leads to an ion fraction that decays with increasing density. For the densest 
(njjj > 10^ cm~^ ) most evolved cores the ultimate expect ation is that Dj^ is 



the dominant ion (|Roberts et al.ll2003l : IWalmslev et"all 120041 '). Observations and 
models now suggest that in at least some objects the ion fraction in the core 
center is quite low, Xe ~ 3 x 10~^, implying that the gas is locally decoupled from 
the field (i.e. the timescale for ambipolar diffusion is compara ble to or less than 
t he free-fall time. jCaselli et al ' '2002c^ : Ivan der Tak et all |2005| ) . 

Walmslev et all tOOi ) and .Flower et al.l ^200^ discussed ambipolar diffusion 



timescales in the context of depletion, providing detailed expressions: 



(5) 



where G is the gravitational constant, i implies a summation of ionic species, 
fiin = mimHi/i'fni + ^-^2) is the reduced mass with H2 the neutral collision 
partner. The rate coefficient for momentum transfer between ions and neutrals 
i s assumed to be < av >in= T^e{a/ Hin)^'^ , with a as the polar izability of H2 
(IOst,erbrocklll96lh . If HCO^ is the dominant ion then Tad ~ 1-7 x lO^^a^g. An 
essential feature is that the timescale is proportional to the square root of the 
reduced mass implying that the ionic composition also influences the diffusion of 



magnetic flux (the timescale is 



lower if 



dominates and 25% for Dg" 



)• 
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These works highlighted the potential of negatively charged grains to alter the 
degree of ionization and also as limiting factors in deuterium fractionation, which 
can be influenced by the changing size distribution as the result of coagulation. 

4.5 Key Tracers and the Possibility of Complete Heavy Element 
Depletion 

Starless cores have revealed a dynamic and evolving chemical structure that high- 
lights key tracers that probe different depths and evolutionary times. One open 
possibility is that all heavy elements are depleted in core centers prior to stellar 
birth. This has important implications for our ability to study the process be- 
cause the heavy element species (e.g. CO, CS, N2H+, ...) have permanent dipole 
moments and can be observed via pure rotational transitions in regions of the at- 
mosphere with good transmission via heterodyne techniques at mm wavelengths. 
If the timescales of star formation are such that heavy element depletion occurs 
prior to the formation of a point source then our ability to study the entire pro- 
cess will be severely hampered, leaving only hydrogen/deuterium molecules with 
dipole moments as probes (H2D"'", D2H"'"). 

At present the evidence for total heavy element depletion is inconclusive and 
focused on the question of whether N2H+ is seen to deplete. This molecule is key 
as it also is the lik ely formation route of NH .cj, and if it deplet es then ammonia 
should follow suit. iBergin et all (|2002 ') and ,Maret et al.1 (|2006l ) find that N2H+ 



does drop in abundance towards the center of the B68 core. In contrast, the N2H+ 
and NH3 emission in LI 544 when ob served in high resolut ion appears to coincide 
with the continuum emission peak ( Williams et al. 19991 . Crapsi et al. 2006, in 



prep.). Thus in the two best studied sources we find contradictory solutions and 
total heavy elemental depletion remains an intriguing possibility. 

Given the gains in our understanding of the chemistry and how it changes 
with core evolution, in Fig. [12] we provide a schematic to summarize our current 
knowledge. For less evolved sources (71^=0 < few x 10^ cm~^) the effects of 
depletion are likely minimized and commonly observed molecules such as CO, 
C^^O, CS, and HCO^, remain viable tracers. For evolved sources {nr=o > few x 
10^ cm~^) species such as N2H''~, NH3, N2D+, and DCO"^ appear to be the best 
molecular probes of core gas, perhaps even to the dense center. Regardless of 
whether total heavy element depletion has occurred, both H2D+ and D2H+ are 
unambiguous tracers of the densest gas. In contrast, emission from C^'^O, CS, 
H2O, and HCO"'" preferentially sample the core edge. 

5 INFRARED DARK CLOUDS 
5.1 Background 

Infrared dark clouds (IRDC) were discovered serendipitously in mid-infrared sur- 
yeys of the ga.lactic plane by the Infrared Space Observatory's ISO GAL survey 
( Perault et al.l[l99^ rl and the Midcourse Space Experiment (MSX; lEgan et al 



11998114 Both surveys revealed a population of objects that were dark at mid 



infrared wavelengths with opacities ranging from 1-4 at 8 and 15//m (I Carey et al. 



■^ISOGAL consisted of a 7 and 15 /im survey of 12 degrees in the plane interior to \l\ = 45°. 
''MSX surveyed from / = 91 — 269° with 6 ± 5° in 5 filters sampling between 4-25 fxra. IR 
dark clouds are best seen in MSX Band A images. 
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19981 : iHennebelle et alJboOll l. The objects appear to be quite filamentary and 
comparison with IRAS suggests that the objects were dark from 7 - 100 fim and 
"represent a p opulation of heretofore unrecognized popu lation of cold, dense, iso- 
lated clouds" te gan et al]ll998l : IHennebelle et alJl200ll ). A compilation of MSX 
dark objects finds over 10,000 sources with kinematic distances between 1-8 pc 
with a distribution that peaks towards prominent s tar forming region s, spiral 
arm tangents, and the 5 kpc Galactic Molecular Ring ( Simon et alJioO^ ). Fig. [13] 
presents one example showing the Gil. 11-0. 12 IRDC as seen by the Spitzer Space 
Telescope and in sub-millimeter continuum emission. 

Much of the interest generated by this discovery has focused on the formation 
of high-mass stars and stellar clusters. This review demonstrates many of the 
gains in our knowledge of the formation of low-mass stars through the isolation of 
starless cores as a unique sample to study the initial conditions of star formation. 
In the case of clustered star formation the nearest region is Orion at ~ 500 pc 
and it is unclear how representative Orion is when compared to star formation 
in the inner galaxy, where most of the molecular mass resides. IRDC's provide 
the capability to examine the early stages of star cluster formation beyond the 
local solar neighborhood. The primary issue for massive star formation is that 
the timescales are much faster when opposed to the formation of solar mass 
stars. Thus there will be fewer pre-stellar massive objects present at any given 
time to isolate. Moreover, since massive stars are born in more distant clouds 
with stellar clusters t he chances of confusion are greater (for a nice review, see 
Garay &: Lizanol 19991 ). The discovery of infrared dark clouds provides a new - 
greatly expanded - sample to search for the precursors to massive stars. Without 
constraining observations, theories of massive star formation are divided between 
models that suggest that high mass stars form by accretion of an extended en- 
velope or competitive accretion of cluster stars for remaining u nbound gas, with 



stars located near the centr al potential obtaining higher masses (jKrumholz et al 



20051 : iBonneh fc BatellioO^ ). A discussion of the theoretical issues is beyond the 
scope of this review and the reader is referred to the cited references and the 
contribution by Yorke &: Zinnecker in this volume. 



5.2 Physical Properties and Implications 

The properties of infrared dark clouds have been estimate d using a variety of 
techniques that have been outlined in §3 andlEvansI (|l999l ). From observations 
of CO, NH3, CH3C2H, and H2CO temp eratures are estirnated to be below 20 
K with densities in excess of 10 cm" -3 (jCarev et al.lll998l : iTevssier et al.ll2002l : 
Pillai et ahlboOfil ). Surveys of mm/submm continuum emission and N2H"'" con- 
strain typical core masses of 100-1000 Mq and sizes of i? ~ 0.3 — 0.7 pc which 
are comparable to those of molecular cores containing ernbedded high-mass proto - 
stellar objects. ( Sridharan et al. 20051 : Ragan et al. 20061 : Rathborne et al. 20061 ). 
The cores are dominated by turbulence with velocity linewidths of ~ 2 — 3 km 
s~^. These linewidths are systematically less than seen towards cores with em- 
bedded protostars, which is similar to th e dependence observed in local clouds 



(jSridharan et al.ll2005 ': 'Raga n et al.ll2006ll. The core mass spectrum has a power 



law slope of 2.1ib0.4 ( Rathborne et al. 20061 ). which is comparable to the Salpeter 
IMF; however, it is likely that these objects will fragment on scales below the 
observed resolution (11")) and this value is also consistent with the CO mass 
spectrum of molecular clouds (§2.3.4). Based on current evidence, it is clear 
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that these sources hkely represent the birth sites of stellar clusters and that a 
sub-sample must be forming massive stars. 

The implications of IRD C's for galactic star formation have been discussed 
by iRathborne et al. who estimate that the star formation rate in these 

objects is ^ 2 M(^/vr. This is comparable to the total Galactic rate of 3-5 
Mo/yr (jPrantzos fc AubertJ [l995l ) and hints that this population is responsible 
for current star formation in the Milky Way. The full promise of these objects 
will be realized as Spitzer galactic plane surveys (GLIMPSE and MIPSGAL) 
are analyzed and sensitive high resolution instruments are brought to bear (e.g. 
CARMA, PdB, eVLA, ALMA) to provide the kind of spatial resolution that is 
currently available towards the Orion Molecular Cloud. Such observations, when 
combined with the techniques developed for local starless clouds described in §3- 
4, will ultimately aid in constraining the fragmentation of molecular clouds and 
the eventual formation of stellar clusters including massive stars. 



6 SUMMARY AND FUTURE PROSPECTUS 

Cold Dark Clouds are a fundamental component of the galactic interstellar medium 
They represent the nearest sample of star-forming molecular material, and as 
such, they play a key role in unraveling the still mysterious process via which 
stars and planetary systems are born. As this review has shown, the study of 
dark clouds has seen some significant progress but still contains significant areas 
of uncertainty. 

On the large parsec-scale, high resolution maps of optical and infrared extinc- 
tion have provided a needed "big picture" look at cloud structure, and have de- 
termined that a significant amount of material lies at low extinction ( Ay < 3 — 5). 
Complementing these data are maps of molecular emission where, aided by im- 
proved sensitivity and mm-wave heterodyne arrays, we can probe deeper using 
^^CO and C^^O than allowed by optically thick ^^CO lines. These data have 
highlighted the filamentary and irregular structure of molecular clouds, sugges- 
tive of a transient nature. The velocity structure, on the other hand, tells a tale 
of the dominance of turbulence, and sophisticated techniques to characterize gas 
motions have constrained the spectrum of turbulence allowing us to move one 
step beyond Larson's Laws. Magneto-hydrodynamic models have matured to 
the point where they can produce realistic turbulence, and numerous simulations 
have shown that the decay of MHD turbulence is rapid and unavoidable with- 
out a continuous source of replenishment. This raises an important challenge to 
established theories of long-lived clouds (> 10 Myr). In addition to turbulence, 
static magnetic fields are a potentially important ingredient of cloud physics. 
Measurements of this field show that in some instances the polarization vectors 
follow cloud structure, but the measurement and interpretation are complicated 
by various effects. 

In regions where turbulence has dissipated, dense self-gravitating cores form. 
By mapping their molecular line and dust continuum emission, we have found 
that cores are characterized by a central region with near constant density and gas 
temperature (to within a factor of 2), together with subsonic motions. Bonnor- 
Ebert spheres often provide accurate fits to the observed emission, and although 
they represent a simplification of the core physics, they should be considered the 
best available starting point for numerical models of core evolution. As in the 
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case of the large-scale, the role played by the magnetic field is still uncertain. 
In addition, detailed radiation transfer analysis suggest the presence of motions 
along the line of sight. These motions may have an inward component that is 
associated with core condensation, but often contain additional contributions of 
an unknown origin. Core-envelope velocity differences, pulsation, or rotation may 
be responsible for them. 

Knowledge of core physical structure has revealed the complex chemical compo- 
sition of starless cores. Observations demonstrate that chemical inhomogeneities 
are an inherent dense cold core characteristic that needs to be considered when 
interpreting molecular line data. The combination of sophisticated radiation 
transfer codes and realistic core physical structures allows, for the first time, 
the investigation of line of sight chemical gradients. As a result, it has been 
found that gas-grain interactions dominate the chemistry of cores, and this has 
a number of consequences. The primary effect is the freeze-out onto grains of 
important gaseous species, such as CO. This removal of CO from the gas phase 
leads to the production of new species, like N2H+ and NH3, in the denser gas 
where these effects are dominant. In addition, the deuterium fractionation of 
the core is greatly enhanced, as multiply-deuterated molecules are produced via 
reactions that begin with H2D"'", 0211+, and D^. Thanks to this new understand- 
ing of the core chemistry, the observed abundance profiles are beginning to be 
successfully reproduced by models that combine chemistry and dynamics of core 
contraction. Using these models we can now isolate molecular probes that target 
specific regions along the line of sight. For example, two classical probes of core 
structure, CO and CS, are now understood to be only tracers of the low density 
cloud surface. Deeper layers of the core are traced by N2H+ and NH3 (along with 
their deuterated daughter products), while the center is unambiguously traced 
by H2D+ and D2H+. These chemical signatures change as a function of core 
evolution and may help to trace the history of gas contraction. 

Within this review we have highlighted areas where progress is needed, and here 
we discuss some areas where breakthroughs could occur using future instruments 
and advance techniques. For the most part we are concentrating on issues that 
can be addressed via observations. 

• There is considerable need to connect the molecular cloud component with 
the atomic diffuse interstellar medium. This may illuminate the question 
of cloud formation, cloud destruction, and the origin of the turbulence in 
molecular clouds, and it will require the combination of mapping observa- 
tions of H I, extinction, [C II], and CO emission. 

• To address the issue of cloud fragmentation and the potential breakdown in 
hierarchical structure, it is required that clouds are mapped with a uniform 
tracer from the large to the small scales. The most likely tracer for this 
work would be the dust continuum emission, if observations retain sensitiv- 
ity to large scale structure, and SCUBA-2 offers the best opportunity in this 
regard. Molecular tracers will provide needed velocity information. Unfor- 
tunately, a wide variety of disparate techniques is currently used to analyze 
the rich information contained in 3-D data cubes. A common description 
of cloud structure is clearly needed both to inter-compare observations of 
different clouds and to test numerical simulations. 

• The transition from turbulent cloud to the quiescent core regime is a crit- 
ical element in the evolution towards star formation. In this transition, 
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turbulence dissipates and the cores acquire the flattened density profiles 
reminiscent of Bonnor-Ebert spheres. Whether the observed subsonic mo- 
tions represent a key part of this transition should be explored. For this, 
an improved understanding of the chemistry should help in the selection of 
a number of molecular tracers to probe different layers. To track turbu- 
lent dissipation, a number of cores spanning a range of evolutionary states 
needs to be isolated and characterized. Millimeter observations with ex- 
isting and future radio telescopes, preferably with large format heterodyne 
receiver arrays (e.g. SEQUOIA on LMT), offer the potential to examine 
these issues. 

• A clear advance results from the ability to combine measurements of the 
dust column spanning a range of wavelengths from ncar-IR absorption to 
submm emission to examine the dust properties, in particular the mass 
opacity coefficient. In parallel to understanding how the gas changes with 
contraction there is an open question as to how dust properties change 
as the mantles grow and grains coagulate. Herschel and SOFIA observa- 
tions should aid to disentangle the coupled effects of dust temperature and 
emissivity on dust emission/absorption. 

• Models and observations have mapped a chemical sequence that can track 
the relative evolutionary status of condensing cores. Observations can be 
used to provide a rough indication of the evolutionary state, but improved 
theoretical models are still needed to provide core studies with an accurate 
chemical clock to test core condensation theories. Further testing of models 
against observables (e.g. integrated emission and line width) also requires 
accurate coUisional rates, molecular rest frequencies, measurements of low 
temperature gas phase reaction rates, and studies of the gas-grain interac- 
tion in the laboratory. 

• It is not clear yet whether the centers of cores are completely depleted in key 
elements (C, O, and N) preceding the onset of gravitational collapse, and 
this has important implications in the ongoing search for star-forming infall. 
Herschel, SOFIA, and especially ALMA will be important instruments to 
clarify this issue. 

• The emerging field of infrared dark clouds holds the promise of unlocking 
the secrets of massive star formation, which is an area with considerable 
theoretical controversy. At present it is not certain how IRDC's fit into our 
picture of Giant Molecular Clouds. Are these sources the dense centers of 
clouds in the Galactic molecular ring, and if so, how do they compare in 
properties (mass spectrum, size, clumping/fragmentation, turbulence) to 
the local population of massive star-forming regions exemplified by Orion? 
Systematic studies with high resolution (< 4") are needed to bring the 
study of these objecst to the level of detail that currently exists for local 
clouds. eVLA and ALMA will be important in this regard. 
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Figure 1: Top Pa riel: Photograp hic image of the Taurus Molecular Cloud taken 
by E. E. Barnard ( Barnard 19191 ) . His notes state that, ^^very few regions of the 
sky are so remarkable as this one. Indeed the photograph is one of the most impor- 
tant of the collection, and hears the strongest proof of the existence of obscuring 
matter in space" . Courtesy of the Observatories of the Carnegie Institution of 
Washington. Bottom Panel: ^^CO J=l-0 integrated emission map of the same 
region obtained using the Five College Radio Astronomy Observatory. Crosses 
mark the location of known protostellar objects and the emission color scale 
ranges from 0.5 - 10 K km s~^. Image kindly provided by PF Goldsmith in 
advance of publication. 
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Figure 2: Maps of IR extinction (contours) and submm emission (greyscale) for 
the p Ophiuchus dark cloud illustrating how the dense gas, responsible for the 
submm continuum emission, represents a small fraction of the total amount (and 
an even smaller fraction of the total volume). Most of the gas in the cloud is in 
the form of low density, low extinction gas. Contours are at = 2, 4, 6, 8, 10, 15, 
and 20 mag with an effective resolution of 5'. Figure from Young et al. ( 20061 ). 
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Figure 3: Figure illustrating the various methods used to determine the column 
and volume density of starless molecular cores. Left Panels: K-band image of 
Barnard 68 and plot of the derived Ak wit h 10^' resolution as a function of 
core radius (taken from Alves. Lada &: Lada ( 200ll ) with Ak vs. radius pro- 
vided by C Roman-Zuhiga. In this method the measured quantity is the H-K 
excess which is related to Ay by the extinction law, parameterized by r^^'^ . 
Ay is correlate d to the hydrogen column from UV line measurements, param- 



eterized by / (|Bohhn et al.l 1 19781 : iRachford et~aD 120021 ). Middle Panels: 1.2 



mm dust continuum emission map and flu x vs. radius for L 1544 (taken from 
Ward-Thompson Motte fc Andrei fl999l and iTafaUa et al.l[2002 ). Hy is the dust 
opacity per unit gas mass, p is the dust density, and m the hydrogen mass (cor- 
rected for He). Rig ht Panels: 7 Atm ISOC AM image and opacity vs. radius for p 
Oph D (taken from Bacmann et al. 2000l ). At 7 pm. the emission from polycyclic 
aromatic hydrocarbons (PAHs) provides a bright background and the dense core 
appears in absorption. In this method the absorbing opacity is related to the 
hydrogen column via the dust absorption cross-section, a\. 
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Figure 4: Left: Molecular hydrogen density, n^/j, gas (T^) and dust (T^;) tem- 
peratures, and visual extinction (Ay) as a function of radius for a marginally 
stable spherical cloud with M = 1 Mq exposed to the standard interstellar 
radiation field (Go = 1). Right: Same as left but for a marginally stable 
cloud (M = 1.6 M0) bounded by an external pressure {pext = 10^ cm— 3 K) 
and illuminated by an enhaii c ed ra diation field (Go = 10). Figure taken from 
Galh. Walmslev Goncah^ ^200i ). 
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Figure 5: Right: Map of dust temperature in TMC-IC derived from 450 and 850 
/im SCUBA maps assuming dust emissivity index (3 = 1.5. Left: Map of dust 
emis sivity index, (3 in TMC-IC assuming a constant temperature of 10 K. Taken 
from Schnee &: Goodman (|2005l l. 
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Figure 6: Sample spectra from the starless cores L1498 and L1517B in the Taurus- 
Auriga dark cloud illustrating the variety of shapes and intensities commonly 
observed. The C^^O spectra trace both dense and low-density gas (due to the 
low dipole moment of the molecule), and often present additional velocity com- 
ponents from unrelated gas that moves at velocities of the order of 1 km (red 
gas in L1498). The optically thick CS spectra are sensitive to the relative mo- 
tions of the core and its envelope thanks to the presence of self-absorptions, that 
commonly indicate inward (L1498) or outward (L1517B) motions at the level of 
0.1 km s~^. The optically thin N2H+ and CH3OH lines trace different regions 
of the core depending on their chemistry, and their nar row linewidth s revea l the 
subsonic nature of the gas internal motions. Data from lTafalla et al. I (|2004l l and 
Tafalla et al.] (|2006l l. 
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Figure 7: A deep opti cal image of the dark globule Barnard 68 (top left; 
Alves. Lada k Ladal[200ll l along with contour maps of integrated intensity from 
molecular emission lines of N2H"'" (contour levels: 0.3-1.8 by 0.3 K km s~^), 
C^^O (0.2-0.7 by 0.1 K km s~^), and 850;um dust continuum emission (10-70 
by 1 mJy beam~ ^ ). Mo lecular data, with an angular resolution of ~ 25" , are 
from iBergin et ah IJ2OO2I I and dust emission (angular resolution of 14.5") from 



Bianchi et all (|2003l ) 



Cold Dark Clouds 



57 



n3 
O 

u 

a 

H 
o 

a 

a 



03 
O 



10" 



Dust Temperature (K) 
10 20 30 



40 



~i — I — I — I — I — r~r-| — I — I — I — I — I — I — I — I — I — I — I — r 



ambi-polar \ 




10' 



10'^ 



n 
o 

m 
< 

o 
o' 

3 

H 



o 

P2 



Logio Density (cm" ) 



Figure 8: Plot of various chemical and dynamical timescales shown as a function 
of the molecular hydrogen volume density. Also shown is the density independent 
CO evaporation timescale as function of the dust temperature (which in the plot is 
also not a function of the density) . Relevant equations for each quantity are pro- 
vided in § l4.3.2l and § I4.4.2[ In this plot the cosmic ray timescale is set by the time 
for any change from an equilibrium state to be reset by cosmic ray powered chem- 
istry. The ambipolar diffusion term assumes = 2 x 10~®(ni^2/10^cm~^)~^'^. 
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Figure 9: Plot of chemical abundances as a function of core depth for selected 
species from a che mical and dynamica l model of a contracting Bonnor-Ebert 
sphere (taken from lAikawa et al. 20051 ). All species are shown at a timescale 
where uh (center) = 3 x 10^ cm~^, excluding CO which also is shown at 
njj (center) = 3 x 10^ cm~^. 




Figure 10: (a) Detec tion of H2D+ emission in L1544 (jCaselli et al.l l2003l : 
van der Tak et al. 20051). (b) 4.4c7 detection of D2H"'" in the starless core 
IRAS16293e (IVastel et ai]l2004l l. 
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Figure 11: Derived electron fraction and major ion abundances in the Barnard 68 
core. Abundances are relative to H nuclei. Zones where the ionization is domi- 
nated by UV and Cosmic Ray are delineated (effectively where C"'' is no longer 
the most abundant ion). In this plot is the main ion, carrying ~ 20% of the 
charge with the remainder contained in more complex ions. Taken from Maret 
&: Bergin (2007), in prep. Some uncertainties in these calculations include the 
gas-grain model (i.e. the assumed binding strengths and desorption mechanisms) 
and the cosmic-ray ionization rate. 
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Major Gas-Phase Tracers in Starless Cores 
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Figure 12: Schematic summary of the major gas-phase probes of starless cores as 
a function of depth and density. The scales provided are approximate and were 
estimated from Barnard 68. In this schematic the temperature is assumed to be 
low (< 15 K) below Ay ~ 2 mag. 
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Figure 13: Gil. 11-0. 12 infrared dark cloud as seen b y (a) Spitzer IRAC b and 4 
and (b) 850 /xm submillimeter continuum emission (I Johnstone et al.ll200S ). In 
the Spitzer image the bright infrared background is provided by PAH emission 
within the 8/im band and the cloud appears in absorption. In contrast, the 
submm continuum emission highlights the dense regions seen in absorption in 
the mid-infrared. 



